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ABSTRACT 


This  report  studies  the  application  of  spectral  domain  analysis 
to  the  problem  of  the  high-frequency  diffraction  of  electromagnetic 
waves.  We  first  summarize  the  basic  formulation  of  the  spectral  approach 
and  then  apply  it  in  a unified  fashion  to  the  Important  problems  of  the 
non-uniform  illumination  of  a half-plane  and  staggered  parallel  planes. 
Special  attempts  are  made  to  analytically  determine  the  fields  at  the 
shadow  boundaries  and  numerically  evaluate  them  in  the  transition  regions. 
Comparison  is  made  with  other  high-frequency  asymptotic  techniques  and 
some  unique  conclusions  are  obtained. 

Furthermore,  it  is  shown  that  the  important  task  of  testing  and 
improving  high-frequency  solutions  can  be  accomplished  by  employing  the 
Fourier  transform  of  the  integral  equation  for  the  surface  current  and 
using  Galerkln's  method  in  the  spectral  domain.  Some  examples  of  this 
approach  are  also  Included. 


1 . INTRODUCTION 
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Ever  since  Keller's  pioneering  work  [1]  on  the  Geometrical  Theory  of 
Diffraction  (GTD) , representing  an  extension  of  Geometrical  optics  (GO), 
there  has  been  a tremendous  amount  of  progress  made  toward  the  development 
of  the  concepts  of  ray  methods  and  their  application  to  high-frequency 
diffraction  problems.  One  of  the  principal  attributes  of  the  ray  techniques 
is  that  they  exploit  the  "local"  nature  of  the  field  solution  at  high 
frequencies.  The  principle  of  the  local  field  states  that  the  high- 
frequency  limit  processes  such  as  reflection  and  diffraction  depend  only  on 
the  local  geometrical  and  electrical  properties  of  the  scatterer  in  the 
immediate  neighborhood  of  the  point  of  reflection  and  diffraction  [2].  This 
concept  is  extremely  useful  as  it  allows  one  to  construct  the  solution  for 
the  scattered  field  from  a complicated  geometry  by  isolating  the  so-called 
scattering  centers,  to  compute  their  individual  contributions  from  the 
knowledge  of  the  scattering  properties  of  associated  canonical  geometries, 
and  finally  to  sum  these  contributions  to  generate  the  total  field. 

Although  the  ray  methods  based  on  Keller's  theory  work  remarkably  well 
for  a very  wide  class  of  high-frequency  diffraction  problems,  there  are 
situations  where  the  need  for  refining  these  solutions  becomes  clearly 
evident.  A few  examples  of  such  situations  are:  (i)  the  observation  point 
at  which  the  scattered  field  is  desired  is  located  at  the  incident  or 
reflection  shadow  boundary  of  a diffracted  edge;  (ii)  a second  diffracting 
edge  is  located  at  the  shadow  boundary  of  the  first  edge;  and  (ill)  near 
end-on  Incidence  of  a single  or  multiple-edged  structure,  e.g.,  a strip 
or  an  open-ended  cylinder.  Various  approaches  based  on  so-called  uniform 
theories  or  equivalent  current  methods  have  been  proposed  and  extensively 
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developed  with  the  purpose  of  circumventing  the  difficulties  in  Keller's 
theory.  Typically,  they  attempt  to  preserve  the  spirit  of  Keller's  ray 
concepts  and  correct  it  postfact  with  some  modifications  of  the  original 
formulas  at  the  "trouble  regions." 

In  this  communication  the  authors  develop  an  alternate  interpretation 
of  the  high-frequency  diffraction  phenomenon  which  is  based  on  the  concepts 
of  "spectral  domain"  [3].  In  this  method  the  scattered  far  fields  are 
viewed  as  the  Fourier  transform  of  the  induced  current  on  the  scatterer  and 
are  thus  associated  with  the  spectrum  of  the  surface  current  distribution. 
Some  significant  advantages  accrue  from  the  use  of  the  spectral  domain 
concepts  for  the  derivation  of  the  solution  of  high-frequency  diffraction 
problems.  First,  we  show  that  the  spectral  concept  yields  results  identical 
to  Keller's  formulas  where  they  apply,  but  the  formulas  based  on  spectral 
concepts  continue  to  be  valid  at  shadow  boundaries,  and  caustic  regions, 
and  for  edge-edge  interactions,  etc.,  where  Keller's  formula  runs  into 
difficulty.  Second,  we  demonstrate  that  the  spectral  formulation  can  be 
used  in  a very  systematic  fashion  to  produce  the  higher  order  asymptotic 
terms  for  some  important  cases.  Finally,  we  show  that  the  Important  task 
of  testing  and  systematically  Improving  high  frequency  solutions,  whether 
derived  by  ray  methods  or  spectral  techniques,  can  be  accomplished  by 
working  with  the  Fourier  transform  of  the  integral  equation  for  the 
surface  current  and  by  using  iteration  or  Galerkin's  method  in  the 
spectral  domain. 

Illustrative  examples  and  comparison  with  other  contemporary  techniques 
for  solving  similar  high-frequency  diffraction  problems  are  liberally  in- 
cluded in  the  paper  to  show  the  ease  of  application  and  the  systematic  nature 
of  the  spectral  concept  and  its  usefulness  for  handling  practical  problems. 


2.  DIFFRACTION  OF  A IM.ANK  V'AVK  BY  A HALF-PLANE 


The  problem  of  plane-wave  diffraction  by  a half-plane  has  been 
analyzed  extensively  in  the  1 iterator'*  since  Sommerfeld's  well-known 
solution  in  1896.  The  reader  may  refer  to  the  standard  texts  of  Noble 
[4l,  Born  and  Wolf  [5|,  Mittra  and  Lee  | 6] , and  others  in  which  compre- 
hensive reviews  are  found.  The  principal  reason  why  the  half-plane 
solution  plays  such  an  important  role  in  diffraction  theory  is  that  it 
forms  an  integral  part  of  the  solution  of  a large  class  of  high-frequency 
diffraction  problems  dealing  with  more  complex  bodies.  In  this  section 
we  re-examine  this  classical  problem  from  a new  angle  in  which  the  solution 
is  constructed  in  the  spectral  domain  after  introducing  the  concept  of 
the  spectral  diffraction  coefficient.  Only  a brief  discussion  of  the 
solution  is  presented  here,  mainly  with  the  objective  of  laying  the 
foundation  for  more  complex  problems  to  be  dealt  with  in  the  following 
sect  ions . 

2 . 1 Basic  Formulation 

The  geometry  of  a perfectly  conducting  half-plane  located  at  y = 0, 

X ^ 0 and  illuminated  by  a plane  wave  is  shown  in  Fig.  la.  The  Cartesian 
coordinates  (x,y,z)  and  the  cylindrical  coordinates  (p,<p,z)  are  erected 
at  the  edge  of  the  half-plane.  Angles  are  defined  positively  counter- 
clockwise with  the  range  [-ti,it].  We  let  the  direction  of  propagation  of 

► i A 

the  incident  plane  wave  be  normal  to  the  edge,  i.e.,  k • z = 0.  This 
assumption  changes  the  vector  nature  of  the  three-dimensional  problem  to 
a two-dimensional  scalar  problem.  Furthermore,  the  problem  may  be 
classified  as  in  the  cases  of  E-wave  (nonzero  field  components  E^,  H^,  H^), 
or  H-wave  (nonzero  field  components  H^,  E^,  E^)  by  simply  letting  the 
incident  field  E-field  or  H-field  be  directed  alternatively  along  the 
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z-axis.  Unless  otherwise  spt'cified,  the  rases  of  E-  and  11-waves  are 
treated  simultaneously,  with  t tu’  help  ot  two  symbols  u and  t such  that 


for 

K-w.i  ves  : 

u = F^  , t = - 1 

(la) 

for 

H-waves : 

u = 11  , t = +1 

(lb) 

The  total  field  u^  may  be  split  into  the  incident  field  and  the 
scattered  field  u^  to  give 


t _ I s 

U = U +11 


(2) 


For  a perfect  electric  conductor  the  total  field  u*"  is  subject  to  the 
boundary  condition  u^  = 0 or  9u*’/9y  = 0 for  E-wave  or  11-wave  cases, 
respectively,  on  the  half-plane.  If  one  defines  the  induced  electric 
current  on  the  half-plane  as 


for  F-wave:  J 

z 


for  H-wave:  J 

X 


9u 

9y 


(3a) 


(3b) 


and  uses  the  time  convention  e , one  can  readily  arrive  at  the  following 

equations  using  Maxwell's  equation 


for  E-wave: 


E^*  = io)  p 
z 0 


rO 


(x' Igpfk |p  - x' X I ) dx' 


for  H-wave: 


s 


u 


_9_ 

9y 


x’xl)  dx’ 


(4a) 

(4b) 


where  k = , p and  r are  the  permeability  and  permittivity  of  the 

medium,  respectively,  and  g^(kp)  ” illQ(kp)/4  (H^  is  the  Hankel  function  of 
first  kind  and  zero  order).  The  objective  is  to  determine  .1  ami  u for 

C2 
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the  half-plane  illuminated  by  an  incident  plane  wave.  This  is  done  by 
using  transform  technique  and  employing  the  results  given  in  [6]. 

2 . 2 Spec t ral  Diffraction  Coefficient  and  Total  Field 
Let  us  define  the  Fourier  transform  pair  as 


li(a) 


,(T> 

u(x)e^°*^  dx  = F[u(x)] 

— oo 


(5a) 


and 


r'-Fi  A 


u(x) 


U(u)e  dx  = F ^[U(a)] 


“+iA 


(5b) 


where  A is  a small  positive  number.  The  incident  plane  wave  may  also  be 
written  as 


i(k  x+k  y) 

X y 

e = e 


i ik*p  _ _ ikpcos  (il-(}i ) 


(6) 


where  k = k cos  Si,  k = k sin  Si  and  0 < Si  < n is  the  incident  angle  shown 
in  Fig.  la.  Transforming  (4)  into  the  spectral  (Fourier)  domain  and 
applying  the  well-known  Wiener-Hopf  construction  [6],  one  arrives  at  the 
following 


for  E-wave:  ~ '^X(k^,ci) 

for  H-wave:  ~ ^X(k^,a) 


(7a) 

(7b) 


where  X(k^,a)  is 


X(k^,u)  = 2 


+ ik^  ^ + ta 


a + k 


(8) 


2 

and  Y » /m  - k such  that  ke  y > 0 and  Im  y < 0,  and  a = k . In  this  work 

~ X ’ 

1/2 

unless  otherwise  stated,  and  (•)  are  defined  with  their  proper  branch 
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cut  slightly  below  the  negative  real  axis.  Using  the  transform  version  of 


(4)  and  then  inctirporat ing  (7),  one  finally  obtains 


U 


s 


X(k^,a) 


(9) 


Furthermore,  one  may  notice  that  the  following  equation  has  been  used  in 
the  construction  of  (9) 


F[gQ(kp) ] 


(10) 


Introducing  the  change  of  variables  x = p cos  (f,  y = p sin  <() , 
k^  = k cos  ft,  = k sin  a = -k  cos  ip  and  y = -ik  sin  ip  into  (9)  and 

substituting  the  result  into  (5b),  one  finally  arrives  at 


s 1 
u = 7— 
4tt 


Ugn  (41) 


X(n,4/)  d4> 


for  ' 


E-wave 


H-waveJ 


(11) 


In  the  preceding  equation  41  is  the  complex  angle  defined  on  the  path  F, 
shown  in  Fig.  lb,  and  x(S7,i|')  is 


X(Q,4^)  = X(k  cos  n,-k  cos  i|;)  = + TXj.(f2,i|^) 


(12) 


where 


Xj(n,4<)  = +CSC  — 

r 


(13) 


We  may  notice  that  Xj(’)  ^md  x^(')  have  the  same  functional  form,  l.e., 

CSC  (•)■  This  definition  of  Xj  ■'^tid  x^.  is  closely  related  to  the  definition 
used  by  Deschamps  in  [7].  Clearly  Xj  and  x^.  are  infinite  at  ij;  = il  and 
4;  = -il,  respectively.  These  two  values  of  w correspond  to  the  incident 
and  reflection  shadow  boundaries  appearing  in  the  GTD  technique.  As  a 
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matter  of  fact,  [irecisely  the  angular  part  of  Keller's  diffrac- 

tion coefficient,  when  w is  replaced  by  the  observation  angle  (f . Although 
X tends  to  infinity  at  the  shadow  boundaries,  it  does  not  mean  that  the 
field  itself  is  also  infinite  as  Keller's  GTD  predicts.  Instead,  the 
correct  value  of  the  field  is  obtained  from  (11),  wiiich  is  always 
bounded.  To  distinguish  it  from  Keller's  coefficient,  which  is  associated 
with  the  diffracted  field,  we  will  refer  to  x(i2,it;)  as  the  Spectral 
Diffraction  Coefficient  for  the  tialf-plane.  This  terminology  is  chosen 
since  xCli.'i')  is  associated  with  the  spectrum,  or  equivalently,  the 
Fourier  transform,  of  the  induced  current  and  appears  only  inside  the 
kernel  of  the  plane  wave  spectrum  representation  for  the  field  and  not 
directly  in  the  form  of  a factor  multiplying  the  incident  field  as  in  the 
case  of  Keller's  representation. 

We  may  further  use  (4)  and  (7)  and  introduce  the  spectral  coefficient 
of  the  physical  optics  field  as  the  Fourier  transform  of  the  physical 
optics  induced  current  to  arrive  at 


for  E-wave:  X*^°(k  ,a) 

X 

for  H-wave:  x'^°(k^,a) 


2k 

iL 

a + k 

X 

2i/g  - k 
a + k 

X 


(14a) 


(14b) 


The  application  of  the  change  of  variables  used  in  (11)  allows  one  to 
express  (14)  as 


where 


x^°(fi,'(/)  * +ctn  . 


(15) 


(16) 
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It  is  worthwhile  to  mention  that  x . defined  in  the  following  equation, 
is  bounded  at  the  shadow  bmindarles 

- x'^^’dl,'!')  • (17) 

could  be  called  t he  f r inge  d i f f rac  t ion  coefficient  and  it  is  used 
for  the  scalar  aperture  diffraction  problems  [8]. 

For  the  problem  at  hand,  i.e.,  incident  plane  wave,  the  spectral 
integral  (11)  can  be  expressed  exactly  in  terms  of  the  Fresnel  integral, 
viz . , 


ikpcos  (S2-<()) 


F(-f,.)  + T 


ikpcos(ii+ii!) 
e FfC  ) 

r 


(18) 


where  the  Fresnel  integral  F is  defined  as 


F(C)  = 


-in /A 

^ i 


. 1 
It  . 
e dt 


(19) 


and  its  properties  are  discussed  in  [9]  and'  [10].  Furthermore,  Cj  and  are 


C.  = +/2kp  sin  - ^ . 


(20) 


L'sing  the  analytic  continuation  argument,  one  can  sliow  that,  for  complex 
angles  of  Incidence,  (18)  is  still  the  proper  solution  of  the  diffraction 
problem.  In  this  context  0 is  replaced  by  the  complex  angle  w which 
follows  the  path  [ (i“,0)U(0, 11)11(11, -i“)  ] , in  the  complex  w-plane  to 

cover  the  infinite  spectrum  of  incidence  angles. 

In  reviewing  the  material  presented  In  this  section,  we  note  that 
its  principal  contribution  has  been  the  introduction  of  the  spectral 
diffraction  coefficient,  which  — In  turn  — Is  shown  to  be  associated  with 
the  integral  representat ion  of  the  scattered  field.  The  equivalence 


between  the  CTD  results  and  those  derived  from  the  spectral  representation 
for  observation  angles  not  close  to  the  shadow  boundaries  can  be  easily 
established  by  substituting  the  asymptotic  expansion  of  Fresnel 
into  (18).  In  the  next  few  sections,  we  will  illustrate  the  broad  nature 
of  ttu-  spectral  concept  and  its  versatility  of  application  by  considering 
more  general  incident  waves  and  complex  structures  than  the  half-plane 
illuminated  by  a plane  wave. 
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3.  DIFFRACTION  OF  AN  ARBITRARY  INCIDENT  FIELD  BY  A HALF-PLANE 

Recently  tliere  has  been  much  interest  in  the  problem  of  predicting 
the  characteristics  of  antennas  mounted  on  complex  structures,  such  as 
aircraft,  and  in  the  accurate  calculation  of  the  radiation  pattern  of 
reflector  antennas  used  for  high-frequency  communication.  Accurate 
solution  of  these  problems  requires  a thorough  understanding  of  the 
scattering  process  of  an  arbitrary  incident  f ield  impinging  on  a diffracting 
edge.  Though  the  problem  of  diffraction  by  a half-plane  due  to  an  incident 
plane  wave  has  been  studied  extensively  in  the  past,  it  has  not  been 
analyzed  in  depth  for  the  case  of  an  arbitrary  incident  field.  The  first 
rigorous  m-athemat ical  treatment  of  the  problem  was  given  by  Carslaw  in 
1899  where  he  considered  the  very  special  case  of  the  incident  field 
due  to  an  isotropic  line  source.  The  latter  problem  was  also  studied 
by  Clemmow  [11],  who  provided  much  insight  into  the  subject  through  the 
use  of  plane  wave  spectrum  analysis.  In  his  paper,  Clemmow  presents  a 
good  review  of  the  existing  literature  up  to  the  year  1950.  A compre- 
hensive study  of  the  half-plane  diffraction  due  to  an  Incident  plane  wave 
and  an  Isotropic  line  source  is  given  in  Chapter  11  of  the  book  by  Born 
and  Wolf  [5].  Khestanov  [12]  derived  an  integral  representation  of  the 
field  for  the  diffraction  by  an  arbitrary  incident  field  and  presented  a 
formal  solution  of  the  problem.  In  his  construction,  no  details  regarding 
the  behavior  of  the  field  at  the  shadow  boundaries  were  provided. 

Recently,  Boersma  and  Lee  [13]  have  employed  the  formalism  of  the  Uniform 
Asymptotic  Theory  (UAT),  which  is  based  on  an  Ansatz  develop'd  by 
Ahluwalla,  Lewis  and  Boersma  [14],  to  obtain  the  asymptotic  solution  of  the 

field  diffracted  by  a half-plane  due  to  a line  source  field  to  the  order 
-3/2 

of  k of  the  Incident  field. 
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Our  objectives  in  this  section  are  the  following:  (i)  to  present  an 

integral  representation  of  the  incident  field  and  study  its  asymptotic 

behavior;  (ii)  to  derive  tlie  radiated  field  behavior  of  a multipole 

line  source;  (iii)  to  use  spectral  domain  analysis  (STD)  in  a manner 

discussed  in  the  previous  section  and  obtain  a general  formulation  of 

the  problem;  (iv)  to  use  asymptotic  techniques  and  construct  higher-order 

-5/2 

asymptotic  expressions  for  the  total  field,  including  the  order  k ; 

(v)  to  employ  numerical  techniques  and  determine  the  field  in  the 
regions  where  asymptotic  techniques  cannot  be  used  conveniently;  and 

(vi)  to  compare  the  results  with  those  of  other  asymptotic  techniques  and 
draw  some  unique  and  useful  conclusions. 

3 . 1 Incident  Field-Asymptotic  Expansion 
3.1.1  General  Case 

As  shown  in  Fig.  2a,  a perfectly  conducting  half-plane  at  x < 0 and 
y = 0 is  illuminated  by  an  arbitrary  incident  field  with  the  following 
spectral  plane-wave  representation 


ilcp^cos(w-lil»^l) 

P(o),k)  e d(D,  -n  < (^^  < TT 


(21) 


where  = 0 denotes  the  phase  center  for  the  incident  field  when  its 
Fourier  transform  is  taken  along  the  dashed  line  shown  in  the  same  figure. 
P(w,k)  is  known  as  the  pattern  function;  it  is  assumed  that  this  function 
is  slowly  varying  in  terms  of  a large  k and  has  the  desired  differentia- 
bility property.  Furthermore,  the  integration  path  in  the  complex 
plane  u>  is  shown  in  Fig.  2b  and  its  characteristic  is  discussed  In  [6]. 

As  a special  case,  if  we  let  P(w,k)  = 1 , we  obtain  from  (21) 


1 1 „(1),.  , 1 f 

* ■ 4 “o  Jr 


f lkp^cos(ur-|4^| ) 


dll) 


(22) 
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which  is  proportional  to  the  field  of  an  isotropic  line  source. 

For  larRe  values  of  k (high  frequency),  one  is  usually  interested  in 
determining  the  asymptotic  expression  of  (21).  To  do  this,  one  first 
assumes  that  P(ui,k)  can  be  expanded  in  terms  of  an  asymptotic  series  in 
the  following  fashion 


P(u,k)  - I 

n>0 


(lk)"“  P^(u) 
n 


(23) 


moreover,  P(w,k)  does  not  possess  any  pole  singularity.  The  path  is 
next  deformed  to  the  steepest  descent  path  SDP  passing  through  the  saddle 
point  u)  = as  shown  in  Fig.  2b.  On  this  path,  the  relation 

Re  [cos  (ti)  - l<f^|)]  = 1 holds;  thus  one  may  introduce  the  following  change 
of  variable 


C 


/2 


(24) 


in  which  t is  a real  variable  taking  the  domain  Substituting  (24) 

into  (21)  and  using  the  fact  that  t is  now  a real  variable,  one  finally 
arrives  at 


u^(Pl.^i) 


lkp^+iw/4 
e 

2/2  IT 


G(t,k) 


(25) 


where 


G(t,k) 


P(w,k) 


sec 


(26a) 


in  which  w is  replaced  by 


J5 


I - i n / 4 


2 Arc  sin 


r 


t + 


, ? 


= + i t'n  ' 1 + it‘  j;  ( 1 + + it'/2  + 

I 

as  t < 0 . 


(2hb) 


Using  (2i),  the  asymptotic  expansion  of  (2ha)  may  lie  written  as 


_ ai  - I (J 

Cft.k)  = I (ik)  ; G^(t)  = P_^(w)  ' 

m=0 


m m 


sec  — — 
m 2 


(27) 


where  (26b)  is  also  implied.  The  complete  asymptotic  expansion  procedure 
[l5l  is  used  for  the  asymptotic  evaluation  of  (25).  In  this  procedure, 


one  first  expands  G (t)  in  a Tavlor  series 
m 


- c^"\o) 

^ / \ r m n 

» ■ Jo  '■<"  * ' 


(28) 


where  C^'^^O)  = -- — G (t)' 
m n m 

at 


and  r is  the  Gamma  function.  Then  one 


t=0 


substitutes  (28)  into  (26a)  and  (25)  and  uses  the  result 


/ t 


jj  j + n)/2]  for  n even 

e ^ dt 

I 

for  n odd 


(29) 


to  finally  arrive  at 


u^(Pl,<>l)  - - 


Ikpj^+iit/ 4 


I (ik)-"  I 

2/27tk  m=0  n»0  /it  T(2n  + 1) 


. G<2n)(0)p--l/2 

n 1 

lkp.+lit/4  . _2 

-*■  2^7^  I (Ik)-"  (O.k) 

n-0 


(30) 
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It  wo  I'omhiiK'  tlio  ordi-r  ot  k ' s in  (iO),  wo  r.iav  oxpross  Iho  incidont 
t io  I d as 


lkpj^  + lTi/4 


IkPj+lw/A 


u (p^,<>^)  - 2/7i¥  ^ iTT^ 


m=0 


Vv’hic'h  has  a rav  fiold  iliaraotor  ist  ii-  [16],  i.o.,  satisfying  the  hunehorg- 

Klino  expansion.  Our  task  is  now  to  determine  X^'s  in  terms  of  P 's. 

m m 

From  (d7)  and  (26b),  the  following  relations  are  easily  constructed: 

G„,(0)  = P (U  I) 
m m I i ' 


.(1), 

J i 

m 


(dtl 

'doj/ 

du 

0) 


- UJ 

PtnCo:)  sec  ^ 


““Uj  I 


G^"^(0)  = d" 
m 


-1 


P (oj)  sec  r 

m z 


ihn 


(32) 


where  D = f— I ~ . Simplifying  the  above  relations,  one  obtains 


r 


G^^^(O)  - -2i 

ID 


1 

^ + I 

, Z 4»  m 
duj 


G^^^(O)  •=  -4 
n 


“■■Uil 


dS^(a))  ^ d\(ui) 

in  . j rn  • ^ r*  / \ 

r~  2 ~~n~  ^ 16 


du) 


du) 


(33a) 


(33b) 


The  higher-order  terms  may  also  be  determined  in  the  same  fasliion. 

Substituting  (31)  into  (30)  and  comparing  the  results  witli  (31),  we 
finally  obtain 


j 


ZjCpj.*,)  - PjCUJjp"''^  + 


Po"<UJ)  + 


\ '■o<l*il>  * i ‘'o<l‘il> 

i P'l'd*,!)  *1  p.dPili 


-3/2 


-3/2 

"i 
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-5/2 


I 7 


n4.i) 


(34b) 


n4c) 


whicli  I'l'mp  1 1‘ t of>  t lu'  asymptotii'  evaluation  of  ti  (e.w.)  up  to  and  inrludin;: 

-5  '2 

till'  order  ot  k ' . VJe  not  i ce  Lha  t tlio  abovi'  proeedure  ran  be  systeirat  ica  1 1 y 

followed  to  determine  any  higher-order  terms  in  the  asymptotic  expansion 
of  the  ineidi'iit  field. 

3.1.2  Mu  1 1 i po 1 e hi ny Source 

In  this  section,  we  study  the  behavior  of  the  radiated  field  of  a 
multipole  line  source  and  employ  the  results  of  the  previous  section  to 
construct  the  asymptotic  expansion  for  the  radiated  field.  A multipole 
line  source,  either  electric  or  magnetic,  may  he  defined  as 


N M 


n-0  m»0 


(35) 


where  (n)  symbolizes  an  nth  order  differentiation  and  I 's  are  constants. 

nm 

The  field  radiated  by  (35)  is  proportional  to 


u (x^.y^) 


J(x;.y;>  { 


k/(x^  - xp^  + (y^  - yp^ 


dx’  dy' 


(36) 


where  the  proportionality  constant  is  iiz'()h  for  the  electric  current  source 
(E-wave)  and  for  the  magru-tic  current  source  (11-wave).  Since 


(-l)"(la)" 


(37a) 


till'  liuiru'i  tr:ii.slorm  ol  ( Ui ) aUm>^  t he  x.-axis  niav  be  expressed  as 


lax . 


u (a.y^)  = j u^(x^,y.)  e ^ dx^  = 


N M , , 

I I k'"-"!  (-la)"|?/c.^  - k^) 


N M 

I I 

n=0  m=0 


for  y ^ < 0 • 08 ) 


/2  2 

- I - k = - i k s i n 

info  (^H),  performing  t iu»  inverse  Fourier  transformation,  and  finally 
comparing  tin*  result  with  (21),  we  arrive  at  the  following  interpretation 
for  the  pattern  tum  tion  of  a multipole  line  source* 


P(u,k)  -=  P(u)  \ ]]  I e (cos  u))''(sgn  (p,  sin  u)'" 

n=0  m=0  ^ 


U9j 


As  exi)ected,  the  patti-m  lunction  of  a multipoli*  line  source  is  Independent 
ol  k.  To  construct  the  asymptotic  I'xpression  of  ( Sb)  , oui'  nay  simply  use 
Ml)  and  substitute  I’, ,(■**)  = F(w),  defined  in  (id),  into  (34).  It  is 
evident  that  a closed-form  expression  for  (3h)  may  he  obtained  by  substi- 
tuting (3'>)  into  (3b)  aiui  writing  tin*  final  result  in  terms  of  the 
higher-order  Ilankel  functions. 

3 . 2 Diffraction  h^  a lia  1 f-P  1 ane 

The  geometry  c)f  a h.alf-plane  illuminated  by  an  arbitrary  inciilent  field 


is  shown  in  Fig. 


2a. 


The  half-pl;me  is  illuminated  by  an  .arbitrarv 
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iiu'iiii’iit  I i I I il  with  its  piano  wave  spoctral  rc[ircson  t a L i f)n  p.ivon  in  (21) 
iml  its  pittiTii  rnm  l ion  in  (24).  Tli.-  total  tiohl  n^  mav  bo  split  into  tho 
iiii  iili'nt  t ii  Kl  u'  ami  soattorod  fiold  ii'’  as  in  (22).  Our  p.oa  1 is  to 
dot  <TiTiiiio  u'  with  t lu‘  application  of  tho  spoctral  diffraction  coo  f I i r i on  t . 
To  this  ond,  tho  incidont  fiold  (21)  mav  first  bo  oxiirossod  as 

u\.  of)  = ; i P(m)  > ^ikmos(w-:) 

^r. 

1 

for  fi  sin  f «in  , (40) 


whoro  s and  Si  aro  defined  in  f’ig.  2a.  It  we  substitute  1’ ( ) = -4-’i4(w  - Q) 

into  (40),  the  plane  wave  (6)  will  be  recovoreil.  Ki-proson  t a t i on  (40)  mav 

bo  interpreted  as  a collection  of  plane  waves  with  spectral  density 

i(4r)  ^P(uj)  exp[iks  cos  (ai  - S2)  ] dw  illuminating  tho  half-plane.  Comparing 

(40)  with  (6)  and  using  (21),  one  may  final Iv  express  the  scattered  field 
s 

u as 


1 


Attempts  were  madt"  to  evaluate  (41)  asymptot ica 1 1 v hv  first  replacing  tho 
r integration  by  its  asymptotic  expressit)n  and  then  calculating  the 
intef'rat ion.  Due  to  the  nonun i f orm  nature  of  this  procedure,  proper 
results  could  not  be  obtained  in  a convenient  manner. 


The  inner  integration  over  T in  (41)  may  bv>  ev.aluated  uniformly  in 


terms  of  the  Fresnel  integral  using  (21)  and  (IH).  Once  this  is  done, 
one  may  then  express  u'"’  as 


I 
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(42;U 


wlu-  re 


11 . ( p , ) = , 

1 4 71 

r 


- ^,iksros(„.-i;)  ^,ik,c-cis  (..+  ;)  ^ 

r,  ^ 

1 


In  lilt’  .ibovf  equation  K(-)  h;is  already  been  defined  in  (I9i  .and  ' s are 
obt.iiiu’d  from  (20)  .ifter  is  replaced  bv  oj,  vi:^.. 


f . , = +i''^ko  .sin 


(43b) 


It  is  evident  that  n‘"’(p,'f)  = u^(p,-i}i),  and  that  (42b)  can  be  interpreted  as 
the  superposition  integral.  Our  task  is  now  to  evaluate  (42b)  asymptot i ca 1 1 v 
tor  ditterent  location.s  of  the  observation  and  source  points  according  to 
the  behavior  of  f ^ and  We  first  deform  i.  to  the  steepest  descent 

patli  SDP  (see  Fig.  2b)  through  the  saddle  point  U = givi>n  by 
Re  (cos  (iii  - '.)  I = 1 .and  Im  [cos  (ui  - .T)  | _ 0.  It  may  be  shown  tli.at  the 
contribution  of  the  closing  segments  of  the  integration  path  at  infinitv 
is  Zero.  Since  there  is  no  pole  in  the  integrand  of  (42b),  tlis  equ.ation 
m.iy  be  written  as 


i 


( 

r 

I 


U?  (l  , f) 
r 


4 71 


+ F(iii) 

■'  SDR 


ikscos  (w-i  ) ikpcos  (ai+ij') 
e e 


F(+t".)  d(.) 


(43) 


liquation  (43)  is  the  starting  point  of  our  discussion  for  the 
asymptotic  evaluation  of  the  scattered  field.  Before  de;iling  witli  tliis 
topic  in  the  ni’xt  few  si’Ctions,  let  us  introduce  some  definitions  thtit 
will  m.ake  our  t-  k convenient.  Observ.ation  directions  Ii  = Si  and  ■t  = - 
are  c.alled  ii  it  (SB  ) and  reflection  (RB)  sh.adow  bcnindaries. 
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ri‘S|UHt  i VI' I V , ;is  shown  in  Kig.  2n  . Regions  in  the  vieinitv  of  SI.  and 
KB  , wluTe  _ and  at  ui  = take  some  given  values  for  different  r's, 
are  called  transition  regions  and  are  designated  liv  T*  and  , respee  t i ve  I v . 

1.2.  1 for  a 1 1 ield_  l^r  (}.  / T*  W f*^ 

In  this  region,  the  Fresnel  function  in  (4'i)  is  substituted  bv  its 
asymptotic  series  from  the  following  ex|>ression 

F(C)  ' 0[-Re(Ce“^’'^^)]  + F(0  for  U|  >>  0 (4ia) 

where 

2 Itt  / a 

^ r(m  + , (44b) 


.and  is  tiie  unit  step  function.  It  is  noticed  that  K(i;)  is  an  odd  function, 
i.e.,  K(-  ) = -K(').  One  may  then  exitress  (43)  as 


u®(p,4i)  = u^(p,((i)  + u^(p,i^) 


(45a) 


where 


r 


- , . Ikscos(uj-n)  lkpcos(uj+<|i) 

+ P(u))  e e 


SDP(il) 


-Re  1+f.j  e 
r 


-lit/4 


du 


and 


u^(p,<') 

r 


1 e^*^P  f 
An  J 


-1? 


. Ikscos(tj-n) 
P(w)  e e 


SDP(n) 


(45h) 


P(^^)  du)  , 


Since  u'^(p,.f)  = u'j(i,2ti  - .(, ) , our  attention  will  on  1 v he  focused  on  the 
asym()totic  evaluation  of  u'! . On  the  SIM’  p.ath,  we  introduce  the  change  of 
v.ir  iab  le 


-ln/4  \ 

*’  rl  1 

r ‘ 

ui  ■ 2 Arcsin 

1 - +l  In 

|l  . it'  ’ 

L 

(1  + i)t/ 

1 + it'/2 

1 ^ } 

+ U 


(4n) 


for  t J 0 
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into  ('t’jc)  to  Imvo  till'  inti'j’rat  ion  porfortnod  on  t lu*  roal  lino.  We  notice 

2 

that  cos  (.  - ) = 1 + it^  and  that  both  fn  [•]  and  are  defined  with 

their  branch  (iits  on  the  nep.ative  real  axis.  Kollowinp,  the  same  steps 
used  in  the  previous  section  to  generate  (30),  we  arrived  at 

U^(p.*)  = g(ks)g(kp)  I (ik)'™  [ (Ik)"" 

^ m-0  n-0 


7) 


where 


lkp+lit/4 


" '27^ 


(48; 


and  G^‘^''^'s  are  obtained  from  (28)  using 
m 


- - - r(m  + 1/2)  (2p)““  CSC^®^^ 

ID  IT 


(^) 


sec  —~z-  ^ P(uj) 


(44) 


such  that,  at  the  same  time,  oi  is  replaced  by  t from  (4b). 
The  diffracted  field  may  be  expressed  as 


u^Cp.i^)  * + ^“^^(0,11)  + u^^(p,(fi)  + 0(k  ^) 


(30) 


in  which  the  r.h.s.  of  (30)  is  obtained  by  using  (47)  and  simplifying  the 
result.  After  performing  all  the  necessary  mar i pu 1 a t ions , i.e.,  using  the 
same  steps  as  in  (32)  and  (33),  we  finally  arrive  at 


g(kp)Xj^g(ks)P(fl) 


(31a) 


g(kp)(ikp)"\Jg(ks)P(n) 

-» v^PtO^  4-  rnc 


g(kp) 


x^P'Cn)  + Xj^P"(n) 


(iks)  ^g(ks) 


(51h) 


I 


2 '3 


ii^^(P,<P)  = ^ g(kp)  (ikp)"^X^g(ks)P(fi) 


+ "I  g(kp)  (ikp)  ^ 


-2x^  + 3x^  P(n)  + 3 cos  — ^ " x^P'(i^) 


+ xp'Hfi) 


(iks)  ^g(ks)  + ^ g(kp) 


I X^PC^J) 


+ cos 


n - i 


2xJ  + 3xJ)p’(fi)  + |xi  + 3xJ)p"(n) 


+ 2 cos  XiP"’(n)  + x^P""(n) 


(iks)  ^g(ks) 


(31c) 


who  re 


r 


7 32  + 'Ji 

+c.sr  — - — 


(32) 


.IS  cIcM'inod  in  (13).  The  <iiffr^u'red  field  u m.iy  then  simply  he  constructed 
hv  iisin^  the  following  relation 


u**(p,4i)  " u^(p,4>)  + Ru^(p,i>)  - u^(p,i|i)  + Ru^(p,2ti  - (>)  . 


(53) 


r 

r 


The  expression  of  U|  is  precisely  KcIIit's  CTI)  solution  [1].  Thus,  the 

leading  term  in  (33)  agrees  exactly  with  the  n-sult  provided  bv  Keller's 

(1  iU  for  the  region  where  such  a solution  is  valid.  In  addition,  our 

spectral  an.ilysis  also  provides  higher-order  terms  whereas  GTD  is 

incapable  of  determining  these  terms.  The  expression  for  u^  agrees 

completely  with  the  one  obtained  using  the*  UA'I'  formulation  |13]. 

Kurt hermore , we  have  constructed  to  show  the  ease  in  using  spectral 

analysis  for  the  systematic  determination  of  higher-order  terms. 

To  complete  the  asymptotic  evaluation  of  v/ , defined  in  (2a),  one 

must  also  determine  u^,  given  in  (43b).  First  notice  that 

r 

s cos  (dJ-Si)  +PCOS  (<i)  + ) = p.  cos  (id  - if).) 

r r 


(34) 


wlufi  its  t r i t'.i  1 i lUtpri't  .It  ion  is  shown  in  lif'.  Ja.  (.  , : j .iro  llio 

i i 

iiioni  i n.i  t os  ot  t iio  iil>sor  v.i  t i on  point  si-on  iron  t tio  souroo  point,  .and 

(.  ;^  ) .ii'o  t ho  ooo  rd  i n.i  t os  of  tho  obsorv.it  ion  point  soon  from  tiio  imap.o 

ot  t Ito  sourci'  point  and  is  mo.nsurod  posit  ivi-ly  o 1 ool'wiso  witii  rant;o 

. Substituting  (54)  into  (Abli),  ono  ran  show  tii.it  tho  SDI’  p.ath 

m.iv  bo  dolormod  to  tho  st  oopost  descent  path  goinp,  through  or  ^ , 

i r* 

roS])oc  t i VO  1 V . This  dolormation  is  v.alid  .is  long  .as  f / T or  ; / i , 
rospoc  t i vi’ I V . ilio  .isvmptotic  ev.iluation  of  (A5b)  follows  t lie  same  sti-ps 
.IS  used  tor  (dl)  .ind  tho  result  takes  the  form  (14)  when  is  repl.iced  bv 
I’.  Wo  thoroforo  .arrive  at  the  following  asymptotic  results 


-u  (p^.1),) 


and 


V 


u®(p,0)  =' 


u"^(p  ,<>  ) 
^ r r 


for  n < ^ 

for  n > ^ 


for  -n  < (|i  ^ 
for  -n  > ♦ ^ 


(53a) 


(55bi 


whore  u is  the  field  radiated  from  the  image  point  of  the  line  source, 
rust  om.a  r i 1 y , u + u^  + Ku^  is  called  tho  geometrical  field  in  the  GTi) 
construction  of  the  field. 

3.2.2  Tot  a 1 f i_ejj.j_  .'i  t the  shadow  bound. ir  i c^s  f)  ^ * W 
Along  these  directions,  the  solutions  given  in  (51)  diverge  to 
infinity,  because  x.  and  take  infinite  values  at  T = and  <f'  = -il, 
respectively.  To  overcome  this  difficulty,  we  go  hack  to  (43)  and  rewrite 
it  as  foil ows 


ik(p,‘(1)  = 

r 


i 

4 n 


+ P(uj)  e 


ik(s+p )cos(w-fi) 


/ 


K »^2kc  sin 

I 


Sl)l> 


— — I dw 


(5b) 
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(Sti)  (Knotis  two  ofiuations  i-or  respond  i ntt  to  ikC,  , ) and  nj^((  rcspec  t i ve  1 v . 

't  is  not  iiod  that  t lio  ovalu.it  ion  of  u^(o,-u)  and  uN(,a>  follows  the 
■>.(mo  [)t  i><a‘du  ros  ilovolopod  in  the*  (ire-vious  sootiotj  and  t ho  i r asyrrififotio 
oxp.insion  .an  ho  ohiainod  t rom  (51).  1 ii  t rodiic  i nj;  t lio  oh.anj'.o  of  vari.ablo 

t»u)  into  (5h).  ono  ro.idily  arrivi's  .at 


u"(p.fn)  = ^ ^lk(s+p)  + U/4 

1 4TT 

" 2 
-k(s+p)t  , 

• e dc 


+ p(t)(l  + /kp  t) 


(57) 


wht'ro  p(t)  = !’(.>)  .it  tor  w is  ropl.iood  by  t from  03).  Wo  then  oxp.and 
2 -1/2 

p(t)(l  + it  /2'  in  torms  of  tho  T.ivlor  sorios  .and  substituto  the  result 

into  ()7).  This  equation  o.an  now  bo  ov.tluated  with  the  holp  of  the 
t o 1 1 owinj;  re  1 .a  t ions 


- r t"  dt 

i ^ao 


(58a) 


l2n  " I r(n  + 1/2)  (s  + p)'""^^^ 


(58b) 


^2n+l  + P)  ^ ^2n-l  2^  T (n  + l/2)k  " ^ (s  + p)~^  s'"  . 

(58e) 

Omit  tiny,  .all  tho  in  t ermed  i .i  t e steps,  one  ran  finally  express  (b7)  .as 
toll ows  [ 10 1 

u*(p,D)  - -u^(p,-n) 


- - -|  g[k(s  + p) ] ( P(n)  + (ik)'^ 


+ (Ik) 


-2 


j p"(Q)  + I p(n)j  (s  + p) 

- -1  J ^ 


-1 


j ^ p"(n)  + P(Q) 


(s  + p) 


-2 


- g(ks)g(kp)  P'(fi)  - (ik)"^g(ks)g(kp) 

® 3s(s  + p)^ 


• (P’(fi)  + P"'(f))]  - g(ks)g(kp)(ik)"^  15ps^  + lOsp^  4 3p^ 

30s^(s  + p)-* 

• [4P’(tJ)  + 5P"'(£5)  + P""'(n))  + 0(k“^^^)  . 


(50) 


j 


2h 
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riu'  lirsl  IiTm  ill  t lu'  r.li.s.  of  Ci^)  ran  lu'  idont  i t uni  with  the  noftativi- 

of  I’lK'  halt  of  t hr  inridriit  lirld  whrn  it  is  romparrd  witli  (i4)  for 

Ij  I.  I’  and  1'^  = 1’,  = 0.  Tlir  total  tirld,  driinrd  in  (2),  may 

thru  Iir  ohtainrd  hv  di' t r rm  i n i nn  n^(,  , ) = iu(,i.2f  - : ) from  (sSa)  .md 

r 1 

(>l).  In  (Jh)  t hr  trrris  nfi  to  and  inrindin)’  t hr  onirr  ot  k do  a>;rrr 

witn  t hr  rrsults  y.iviai  in  [lil  whirli  usrd  t iir  lh\i  trrhnii|ur.  It  is  noire! 
that  till'  dominant  asvpi|itotir  trrm  of  t Iir  inridiiit  field  is  of  the  onirr  k 

1.2.  1 Nnmrrii'.il  i val^ition 

In  t hr  transition  rrjiions  T*  and  f'",  ihr  prorrdnros  that  were- 
I'mplovrd  in  t hr  pri'vions  srrtions  for  deriving  the  .isvtnpt  ot  ic  r.xprrssion  of 
t lu’  field  no  l(Mii;rr  applv.  In  this  srit  ion  wr  disriiss  the  niinirriral 
rvaliiation  ot  the  i-xp  fi'ss  i ons  for  t lu  field,  vi/..  (21)  ,ind  (42h),  for 

observation  anplrs  that  I ii'  in  these  angular  regions.  Siiur  the  mirirriral 
lirorrdiire  is  not  rrstrieted  to  ,i  spn  itii  ani'iilar  rr^ion,  wi-  have  rPiplovrd 
it  to  rv.ilu.itr  the  (lit  I fa.  tril  tirld  for  a w id»  r r ni(’,i‘  ot  observation  .iliKlrs 
t h.in  just  the  transition  rrjions.  7hoin>b  these  results  are  not  orderni 
asvm[it  ot  i ra  1 1 V in  terms  (<f  powars  iit  f,  t luv  <li  providi-  .a  pood  basis  tor 
ehrikinp  the  .is  rmpt  ot  i r rt  ■ ai  I t s d.riviil  railief. 

In  their  present  toim,  (21)  .nid  (tjb)  ari’  not  suitable  for  niimrrieal 
intr>>ration  (ju.idr  at  urrs , bn  atisr  ot  tluir  hiphiv  osi  ill.itoiv  kernels.  Ve 
deform  their  paths  to  the  steepest  desenit  paths  and  introdiuc  the  (.iianpe 
of  variables  (2hb)  and  (A<i),  respeel  Ive  ly , in  (21)  .ind  (42b)  to  arrive  at 


lkp^+ln/4 

■ 

^ ^ 2/2T, 


r«D 

P(ui) 

+ P(u.) 

0 

t 

-L 

2 -1/2 

(1  + it  /2)  dt 


Uida) 


where  111  is  replaced  by  t from  (2bb)  and 
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u®(p  ,1^) 
r 


lk(s+p )+in/4 
e 

2/2n 


2 

2 

.00 

r / \ 

[ \ 

0 

+P(w)  e F 

+ P(w)  e F 

k 

\ r'  1 

t ' 

(1  + itV2) 


■1/2  -kst 

e 


dt 


(60h ) 


wl\iTo  . is  ri'pl.'Kfd  by  t from  (4(>).  A sutt'icivnl  romlilion  for  t lu>  iiUo- 
gr.iiul  of  (bOb)  to  bo  boundod  is  tii.it  s p or)S  Cl.  + f ) 0,  respoo  t i ve  1 v , 

for  o'?  or  u'\  Oiio  to  t bo  t ao  t that  tbo  intoprands  of  (bOa)  ami  (bOb) 
havo  smooth  bohaviors  and  dooroaso  very  rapidly  for  largo  values  of  t, 
oonvontional  mimorical  intograt ion  toohniquos,  o.g.,  tlauss  quadratures, 
m.iy  be  omph'yod  to  evaluate  the  integrals  with  good  aecuracv.  Some 
tvpioal  res\ilts  obtained  in  this  manner  are  shown  in  the  next  section 
when-  I Dmp.i  r i sons  are  made  with  ot!u>r  techniques, 
i.  5 Comparison  with  Otlu‘r_  Hi  gh- Frequency  Te^hn  icj^uej^ 

Fver  since  the  classic  paper  by  Keller  [1]  on  the  Ceometrical  Theorv 
of  Uiffraction  (flTD),  numerous  attempts  have  been  made  toward  extending  the 
domain  of  Cil)  to  regions  when*  it  predicts  an  infinite  field  and  is  th<>re- 
fore  invalid.  Most  of  the  forts  in  this  direction  have  been  cc>ncent ra ted 
on  constructing  uniform  formulations,  such  as  UTD,  MSI)  .ind  HAT,  that  over- 
come the  difficulties  i>f  CTI)  at  the  incident  and  reflection  shadow 
boundaries,  and  .at  the  tr.ansition  regions.  The  Uniform  Theorv  of  Diffraction 
(I'TD)  .and  its  Modified  Slope  Diffraction  (MSD)  version  have  been  introduced 
by  Kouyoumjian,  I’.ath.ak  .and  Hwang  1 1 7 ],  [ 18  ] , and  the  Uniform  Asvmptotic 
i'heorv  (UAT)  h.as  been  developed  hy  Ahluw.ali.a,  Lewis  and  Boersm.a  |14|.  The 
l.atter  has  been  employed  extensively  by  Hoersma,  Lee,  Deschamps  .and  Wolfe 
[13],  [19],  [16],  [20]  to  Investigate  a number  of  problems  of  this  type. 

UTD,  MSD  .and  UAT  are  asvmptotic  techniques  based  on  individu.al  hypotheses 
or  Ansatz.  Typically,  there  is  no  general  proof  av.ailable  for  the  validitv 
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or  t lu‘  lomp  U'tfiu’ss  ot  tlu'Ko  Ans.itz  ;iiul  ono  l).is  to  ;ipplv  thorn  to  cortait! 
tost  casos  in  order  to  ost.ililish  their  accuracv.  The  1 Ormu  1 .a  t i ons  ot 
t hi*si’  asvmi'totio  theories  are  reduced  to  t hi'  classical  Somme  r t e 1 d ' s 
ri'siilt  when  the  plane  wave  illumination  on  a half-|>Iane  is  considered. 

In  tact,  Sorami'rfe Id ' s formulation  is  the  basic  foundation  of  all  the 
different  Ansatz  that  have  been  proposeil  to  date.  Kor  more  complex 
situations,  the  validitv  of  the  various  asymptotic  theories  is  typically 
clunked  against  numerical  results  or  experimental  data,  since  analytical 
results  are  often  unavailable  in  closed  form.  Even  for  the  problem  of 
diffraction  of  a half-plane  illuminated  by  a line  source  possessing 
isotri'pic  or  nonisotropic  patterns,  there  is  no  substantial  clieck 
available  for  estah  1 i s!i  ing  the  accur.tcy  of  the  various  asymptotic 
theorit's.  Spec  i f ica  1 1 y , at  the  transition  regions,  it  is  not  known  how 
well  thesi’  asymptotic  theories  compare  with  each  other  or  with  the  exact 
so  1 lit  ion  . 

The  object  of  this  section  is  to  compare  the  aforementioned  asymptotic 
theories  with  the  exact  solution  that  has  been  constructed  using  the 
spectral  analysis  presentiui  in  this  paper.  The  com|iarison  has  been 
carried  out,  both  anal;,  I'cally  and  numerically,  tor  a wide  range  of 
observation  angles  including  the  shadow  boundaries,  transition  regions,  and 
at  the  angles  away  from  these  regions.  For  completeness,  different  source 
locations  and  pattern  functions  have  been  investigated  and  extensive 
results  have  been  derived.  In  the  following,  we  begin  hv  reviewing 
briefly  the  nature  of  the  formulation  of  each  ot  the  uniform  theories, 
and  follow  this  with  the  [iresentat ion  of  the  results  and  conclusions. 

3.3.1  fkirmu  I ation  of  d i f I erent  theories 

All  the  uniform  asymptotic  theories  essentially  approximate  (22)  and 
(42b)  in  an  asymptotic  sense.  It  is  Interesting  to  note  that  as  yet  there 
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is  no  |>roi>t  i v.i  i I ah  I c tliat  justilii-s  t iu‘  use  nt  t lifsc  forms  for  arbitrary 
and  oomi'lox  situations.  In  all  of  t heso  t (-•(.lin  i qius  t lii'  followinj’,  is  true 
fc'r  till'  total  fii'ld 

u ' ( . , >. ) = V (i  , ; ) + t V { , 2 n - ; ) (61) 

wlu’ro  v(i  takes  d_iftertait  forms  for  different  theories.  For  convenience 
ot  in t e rpre t at  ion , all  ol  the  asymptotic  theories  are  schematically  repre- 
sented in  Fig.  i . The  notations  in  this  figure  are  explained  below. 

We  first  give  the  expression  of  v(p,4.)  for  the  GTl)  theory,  which  reads 
as  follows 


v(p,i^)  •=  u^(p,!>)  + u'*(p,i|i)  + 0(k 


(62a) 


where 


u*^(p.^)  - e 


sin  - (n  - 4>) 


and 


■ -CSC  ^ ~ g(kp)u^(s,W) 


(f.2b) 


(b2c) 


In  the  preceding  equations,  0 is  the  unit  step  function  and 

“ g(kpj)l'(t.)  . (b3) 

that  is,  the  first  asymptotic  term  of  (31).  It  is  clear  that  the  CTU 
formulation  fails  at  the  shadow  boundaries,  i.e.,  f = *fl,  and  their 
vicinities. 

The  LTD  formulatiot)  overcomes  this  failure  of  GTl)  at  the  shadow 
boundaries  by  expressing  v(p,(>)  from  (61)  as 

v(p,4i)  ' u^(p,(|i)  + u^^(p,(}i)  + 0(k  (64a) 

R 1)1 

where  u^  is  already  defined  In  (62b)  and  u is 
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u^\p,$)  = 2/^  sgn 


. 4>  - n 

sin  — - — 


/2kps  ..  . -i(n/4+4  ) i 

— R(kp)  c F(4>u  (s,n) 


(hAb) 


In  Ct'-ih)  !'  i .s  t lii'  Krf.sncl  inti'>;r.il  (2b)  ;iiul 


sin 


(()  - n 


/I 


2kps 
p + s 


( bUc  ) 


Vvliiili  niH's  to  zoro  at  t ho  .shallow  houmiarics.  'I'ho  UTl)  f ornui  1 .i  L i on  i ,s  ot  ton 
roiorrod  to  in  tho  litor.itiiro  as  ilio  Kouvoiini  j Ian-Pat  hak  (KP)  formulation. 

Sinco  I'Tl)  is  onlv  suitahio  for  piano  wave  or  isotrojiio  lino  source 
t vpos  ot  inciilont  liohl,  rooontly  a mod  I f I od  version  of  L'TD  that  tmploys 
tho  notion  ot  tho  slope  diffraction  coofficiont  has  boon  introduced  to 
^onorilizo  this  procedure.  This  f ormii  1 ,i  t i on , whicli  is  referred  to  here  as 
MSI),  modifies  (b4a)  as 


v(p,<(>)  •=  u^(p,i}>)  + u^^(p,({))  + 0(k~^^^) 


(65a) 


where 


u^^(p,<C)  •=  u^^P,')’)  + \2/ir  sgn 


sin 


2ksp  p(s  + p) 


■3/2 


sin  (<i  - n)  e ^ F(4)  - 2p(p  -f  s)  ^ cos  ^ — 


g(kp) 


3u 

3^ 


(65b) 


B,n 


In  summary,  UTI)  and  its  modification  MSI)  have  boon  introduced  to  modify  u 
in  tho  (ITI)  formulation  in  such  a manner  .as  to  cancel  the  discontinuity  of  u^. 

Finally,  tho  Uniform  Asymptotic  Theory  (PAT)  has  bi'on  invented  to 
circumvent  tlio  sliadow  l>ound;iry  difficulties  of  the  UTI)  by  writing  v*^  as 


v^(p,4i)  • u^(p,li)  u'^(p,if)  + 0(k 


(bba) 


.F"  'M 
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whoro  u is  llu'  saiiu'  tfrm  used  in  tlu'  (ITD  t’orniu  1 .1 1 i on , i.c.,  (fjlc).  ,;) 

has  till'  torm  tliat  oxaitly  canct'ls  tho  inlinitc  va  1 uf  of  at  the  shadow 
hoiimiar  ii’S.  In  tliis  theory  u*’  takes  the  t Orni 


u‘'(t  ,'t>)  = U-(n  - ) 

1 L 

where  K is  the  lirst  term  of  (44b) 

. . 2 

in/4+i' 

I ' ( • ) = ‘-i 

2/tir. 

ami  C.  wliieh  is  l alli'd  t hi'  detour  funetion,  is 


( 66b ) 


(6he  ) 


4 = SKii  sin  ’ ^ ,'k  j (,  + s)  - I . I ; 

Noll'  that  at  the  shadow  boundaries  the  di'toiir  4 is  xero.  The  higher-order 

-2 

asvmptotie  expression  ol  I'AT  to  the  order  of  k “ (dominant  nsvmiuotic  term 

i . -1/2 

ot  u is  k ) tor  the  probli-m  at  hand  is  given  bv  Boersma  and  I.ee  |li]. 
it  is  worth  mentioning  that  in  all  of  the  above  asvmptotie  teehniques  t he 
field  imident  on  the  diflrartitig  edge  must  be  .1  rav  field  (lo-al  plane 
w.ive).  Somi  difficulties  related  to  the  applieation  ot  in>  idi  at  lields  that 
Jo  not  fulfill  this  criterion  are  discussed  in  Section  *. 

5.5.2  N timer  ica  1 results  ami  discussions 

In  the  previous  sections,  some  results  based  on  analvtii al  and  numeri- 
cal .isymptotie  evaluations  of  (21)  and  (42h)  were  presented  tor  three 
situ.it  ions  point  ot  observation  exactly  .at  the  sh.ulow  bound. tries;  inside 
the  tr.insition  regions;  .ind  outside  the  tratisition  regions.  Addition.il 
results  for  this  problem  of  d i f f r.ic  t i on  of  a nonisotropic  line  sourei-  h.ive 
been  given  in  [1J|.  A comp.irative  study  of  the  various  solutions  .illows 
one  to  deduce  certain  conclusions  discussed  below.  First,  spectr.il 


. j = ►'4“  + p*"  + 2so  cos  (fi  - ) 

( 66d ) 


nacjiBW 


t I'fmii  1 .U  i on  (S'l'i))  i .m  hr  ir  rd  in  .1  vi-rv  svsfin.it  ii  in:l  st  t .1  i I'ln  1 rw.i  rd 
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t'ashirn  in  drfrrminr  .iiiv  niimhrf  nf  li  i f*‘ 'rd>- r t rrm--  in  f hr  .isvmi'ti't  ii 
r xp.ins  i nn  rt  t tlr  li't.il  lirld.  In  ! .11  t , hv  iisiii;’,  tni  . .1;  prii.uli  Wr  li.ivr 
i-rniiint  ril  t hi'  isvmp  t nt  i i rxprrssimi  ot  t tir  [rt  il  tirld  up  [r  1 hr  rrdrr  rt 

k . Srirnd,  I'A'I'  f r rnu  1 .1 1 i on , in  it-.  p,rnrr.ii  vrrsirii,  > in  ilsr  turnisli 

t iir  liiKlir  r-ordi-r  .isviuptotir  rxi'rrssion  ()t  llir  lotil  lirld.  llrwrvrr,  in 
this  nirtlu'd  t hr  hiKhrr-ordri  trrms  in  thr  diltr.ii  trd  field  n*'  nm'~l  hr 
peilriMted  tiirrip’h  .m  rxplirit  ;il>p  1 i r .1 1 i rn  of  the  rd);r  rrnd  i t inn.  Rrsnlts 

lip  tn  the  order  of  k “ have  been  reported  in  | I 3 | , and  t itrv  agree  prrfntlv 

with  the  STI)  results.  Third,  l''ITi  provides  the  i-orrert  .isvmptotir  exiires- 
sion  of  the  field  up  to  the  order  ol  k \ hnt  on  I v tor  .an  isotropic  line 
source.  Furthermore,  I'TD  formulation  does  not  provide  h i giie  r-order 
asymptotic  terms  in  a direif  fashion.  Finallv,  MSI),  which  is  the  modified 
Version  ot  ITD,  .illows  one  to  determine  t lie  asvmptotic  values  of  the  field 
for  tile  iinisotropic  line  source  up  to  the  order  of  k i\gain,  this 

tormul.it  ion  does  not  provide  higher-order  terms  in  a svstematic  manner. 

All  of  the  above  observations  have  been  made  on  the  basis  of  .1 
comp.ir.it  i ve  stmiv  of  the  v.irious  re.sults  directly  at  the  .shadow  hound. tries 
f'  = r = 0)  .ind  far  away  from  the  transitiiui  regions  (f.  • I),'  0)  wiiere 

the  tield  can  he  expressed  in  .1  closed  form.  For  the  s.ike  of  comp  1 1 eness 

we  li.'ive  also  i n vest  i g;i  t ed  the  behavior  of  the  fields  predicted  by  the 
various  theories  in  the  neighborhood  of  the  fr.insition  regions.  For  fhis 
[lurpose  the  field  inti'grals  in  STD,  i.e.,  Fips.  (hDa)  .tnd  (hOb),  have  been 
evaluated  using  an  efficient  numerical  algorithm  discus.sed  in  Section  3.2.1. 
Note  that  the  computed  results  obtaineil  in  tliis  manner  cont.ain  all  of  the 
higher-order  terms.  Numerical  results  have  also  been  obtained,  up  to  the 
r.'inge  of  higher-order  terms  ;ivi  i I ah  I e , using  (ITD,  FTl) , MSP  and  UAT 
f ormu  1 .1 1 i ons  , .as  giveii  in  (h2.i),  (h4a),  (baa)  and  (hha),  respec  t i ve  I v . 


i4 

III  orcior  to  grl  iH'tti'r  insight,  we  liavi'  compared  both  the  total  field 

and  the  scatti'red  field  u^,  deiined  in  (2).  Some  typical  results  are  shown 

and  disiiissi-d  in  I lie  following. 

I'igiiri'  4 sliows  the  total  lield  diffractt'd  by  a liait-plane  i 1 1 umi  iiaLed 
bv  a non  i sot  rop  ic  line  source'  (b-wavi')  with  the  following  pattern  function 
obt  a ined  f roni  ( 19)  . 

I’ ( . ) = cos  '0  + sgn  ( ' . ) s i n . . (fi7) 

It  should  be  noted  tliat  when  (67)  is  used  in  (tiOb),  <f  . is  repl.iced  by 
and  hence  sgn  (ie)  = +1.  Results  are  shown  for  different  techniques  and 
till'  field  was  sampled  in  a 40°  angular  region  about  the  incident  siiadow 
binindarv  (Sb').  The  SB*^  direction  is  shown  in  the  figure,  and  as 
expected,  I'TI)  formulation  fails  to  provide  the  correct  result  due  to  the 
lack  of  pro|ier  slope  discontinuity  compensation.  Since  t lie  incident  field 
is  well-behaved  in  the  entire  rt'gion,  we  have  also  compared  the  scattered 
field  iT'’  for  different  tecliniques  in  Fig.  ' > . Exactly  at  the  shadow 
boundaries,  we  can  employ  our  li  i giie  r-order  asymptotic  results,  derived  in 
C)9)  and  (51),  to  determine  the  scattcreil  field.  i'lie  results  of  these 
calculations  have  been  listed  in  Table  I for  different  orders  of  k.  In  thi' 
s.ime  talile  the  result  of  evaluation  of  (bOli)  lias  also  been  given.  Though 
this  evaluation  implicitly  contains  all  of  the  higher-order  terms  of  k, 
the  agrei-ment  with  the  results  given  in  tlie  table  for  terms  beyond  k ' is 
remarkable.  This  is  iiuleed  a strong  indication  of  the  usefulness  of 
risymptotic  techniques  for  convenient  computation  of  the  diffracted  fii'ld. 

We  have  also  plotti’d  the  results  tabulated  in  the  third  column  of  Table  1 
in  Fig.  5 with  a view  to  exhibiting  thi'  behavior  of  the  asymptotic  solution. 
The  number  of  terms  used  In  the  asymptotic  evaluation  is  indicated  bv 
numerals  inserted  alongside  the  points  representing  the  evaluation  ol  the 
diffracted  field  at  the  shtulow  boundary. 
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Fi^uri-  (i  ilispl.ivs  ilif  bi'li.ivior  of  llio  sr, it  to  reel  field  c.i  I cu  I ;i  t ed  using 
the'  .it’e'reme'ut  iv'U^'d  tei'hniques  in  the  trnnsition  re'gion  iif  a \ialf-plane 
illuninated  hv  an  isotropii  line  source  (i.-wave).  As  expected,  IITI)  and  MSD 
arc  the  same  and  all  the  tis'hniques  are  in  very  good  .agreement.  For 
comp  let  iuiess , the  .asymptotic  expressions  have  been  t.ahul.ated  in  Table  1, 

.and  .are  shoiyn  in  Fig.  h.  Figiin  7 shows  the  pli't  of  the  sf'.attered  field 
from  .1  h.alf-pl.iue  illumin.ited  by  .an  isotropic  line  sourii-  (F-w.ave)  eoplan.ar 
with  t tie  ha  1 f-p  1 .me . For  this  geometry  t tu'  incident  .and  reflected  siiadow 
hound. iries  coincide  with  the  half-plane  itself.  The  scattered  field  was 
s.impli’d  in  .an  angul.ar  regit'n  -10“  < < 10°  .and  nsults  tor  different 

techniques  including  OTl)  .ire  shown.  For  the  sake  of  completeness  the 
.isymptc't  ii'  results  computed  from  (f2.a)  are  also  shown  in  Fig.  for  the 
observ.at  ion  angle  p = 0,  .and  the  oscill.atorv  n.ature  of  this  .asymptotic 
solution  is  evident  from  the  figure. 

To  further  illustrate  t lu’  comp.i  r.i  t i ve  nature  of  the  results  derived 
from  different  theories,  we  employ  t lu-  f ormu  1 .a  t i ons  of  TAT,  I'Tl),  .and  MSD  to 
oht.ain  some  t hree-d  i minis  1 on.a  1 representations  of  the  tot.al  field. 

Figure  8,  containing  the  geomi'trv  of  the  problem,  is  the  three-dimensional 
plot  of  the  total  field  1 1/ | difir.acted  from  a half-pl.ane  illuminated  by 
an  isotropic  line  source  (F.-w.ave).  lo  oht.ain  the  total  field  t/  , we 
employ  the  formulations  of  HAT  and  I'TI),  and  find  that  essentially  identic.al 
plots  are  oht.ained,  .ilthough  numiTical  v.a  I ues  of  t hi'se  curves  differ 
very  slightly  from  each  other.  For  instance,  at  the  observation  point 
X = y = .25>,  UAT  .and  UlT)  formulations  give  |u^|  = .05913  .and  [u^l  = .05918, 
resjiect  i ve  1 y . Figure  9 is  the  counterp.art  of  Fig,.  II  for  the  H-wave  casi>. 

As  expected,  the  H-field  is  discontinuous  at  the  h.a  1 f-p  I .ane . For 
.additional  results  and  three-dimensional  plots  the  reader  is  referred  to 
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Scattered  field  u from  a half-plane  illumlnjted  bv  an  isotropl 
(K-wave)  with  |u^(l,'/2)|  = .0796.  Points  1-6  are  the  r.urber  o 
asymptotic  expansion  of  p.iven  In  Coliirnn  of  Table  1. 
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Till'  mill  t iple-i‘(lgf  d i t 1 r.ir  I i on  prohlfm  irisos  in  ,i  nnmlu-r  ot  |>rai- 1 icn  1 
s i t uat  i iins  and  is  i.  ons  i de  rail  1 v murt'  coni)  1 ii-atod  lliaii  its  sin^lo-edge 
I'ount  orpart  . I lio  comp  1 t-x  i t y ot  t lii  s prohlom  stems  I’rom  liu'  fact  tlial  t he 
lield,  wliicli  im[iin>;i'S  on  tlie  second  islpt'  attir  diffracting  I rum  the  first, 
iMii  no  longiT  be  treated,  in  general,  as  .a  locally  uniform  plane  wave. 

It  has  been  risogni/ed  that  lor  tlu'  problem  of  diffraction  by  a nonplane-wave 
soiircL',  e.g.,  one  with  a pattern  tiinction  as  in  Section  i,  the  result  for 
the  si-atteretl  1 ield  at  the  shadow  boundaries  is  not  just  a simple  moditi- 
cation  of  the  [ilaiu'-wave  solution  multiplieil  by  the  pattern  function  but 
involves  till'  derivative  ot  the  iiattern  function  as  well.  One  might  assume 
tli.it  the  problem  of  mu  1 1 i ]i  1 e-ed ge  ditir.iction  c.an  he  solved  by  a simple 
extension  of  the  p.ittern  function  illumination  of  a single  edge.  That 
this  is  not  the  r.ise,  howevir,  li.is  lu'en  pointed  out  by  dunes  (21  1,  who 
has  employetl  the  l\  i i*ne  r-Hop  1 teclinitpie  to  constnnt  t fie  rigorous  liigh- 
frequeni'V  I'Xp.ansion  ot  the  I'X.ict  solution  tor  the  dittracted  field  by 
sl.agg.i'red  p.ir.illel  pl.ites.  Although  lones  |2l|  did  not  explicitly  obtain 
the  results  for  the  field  .it  the  sh.idov  tiound.i  r i es , hi-  did  examine  the 
uniform  expression  he  derived  in  sotiM'  det.iil  .ind  concluded  that  the 
I'nifonn  Asymiitot  ic  Theory  (FA'f)  could  not  be  directly  (meclian  ica  1 1 y ) 
applied  at  the  shadow  boundaries.  lee  .ind  lloersm.i  |19|  .ind  lioersma  [22| 
have  carried  out  r.atiier  el.ihorate  i n vest  i g.i  t i ons  ot  this  problem  and  have 
shown  that  the  I'AT  can  still  be  used  tor  this  situation. 

In  this  section,  the  (iroblem  <>f  d i f f r.n  t i on  of  a plane  wave  incident 
on  two  staggered  parallel  h.a  1 f-(>lanes  (mu  1 1 i p I I'-i-dge ) is  Investigated  ;ind 
an  analysis  is  presented  which  has  the  tol lowing  featur»-s.  First,  a 
general  representation  of  the  field  alter  successive  diffraction  by  two 
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ha  1 1 -p  liini's  is  ;;ivt'n  in  ti-rms  i)t  a ih'iihlf  cnmpli’X  intfgral.  The  analvsis 
is  larrii’il  out  in  t ho  Fourior  truisti'rm  domtin  and  is  h.isod  (>n  t lie 
applioition  ot  t ho  rosiilts  i;ivon  in  So<  tions  J and  i.  Soomid  , asymptotio 
t (.'ohii  i(|uos  aro  omplovod  to  di'tormino  t ho  total  tiold  to  t ho  order  of  l< 
in  a rip,orous,  vot  s t ra  i «h  t t orwa  rd  fashion  for  some  o.asos  of  jirii  tical 
intorost.  Those  spoiiil  i asos  are  (see  Fig.  10):  (a)  t lie  edge  of  the 

t ('p  half-piano  lies  at  or  awav  from  the  shadow  houndarv  of  the  bottom 
half-piano,  i.e..  the  diroitlv  illuminated  one;  (h)  the  observation  point 
lies  at  or  away  from  the  shadow  boundary  i>f  the  top  half-plane.  In  all 
of  these  oases  the  final  results  are  expressed  in  a compact  and  useful 
form. 

4 . 1 Ba  s i 0 F o rmu  1 ;i^t  i i2n 

The  geometry  of  a pair  of  staggered,  perfectly  londucting,  parallel 
half-planes  illuminated  hv  a plane  wave  is  shown  in  Fig.  11.  The 
bottom  and  top  planes  are  labeled  1 and  2,  respectively.  The  edges  are 
sc’parated  by  distance  s,  and  the  angle  between  the  line  joining  the 
edges  .and  the  x-axis  is  . I.et  the  incident  plane  wave  with  unit 
.amplitude  and  incident  .angle  . take  the  following  form  at  point 

ikD,cos(Q,-<t.  ) 

u ' (p  I , if.  I ) = e . (b8  ) 

In  this  work,  wi-  neglect  the  interaction  between  planes  1 and  2,  because 
it  is  assumed  that  s/X  is  a large  number.  The  total  field  diffracted  by 
plane  1,  u ^ , may  now  he  expressed  as 

U|  = u'  + u^  (b9) 

where  u”  is  the  scattered  field  due  to  the  induced  current  in  plane  1. 

From  (11),  the  scattered  field  u^  can  be  written  as 


OBSERVATION 


It  :v<  i>o^>  i t ions  of  tl)e  oJyo  of  Ltu  top  h.  If-plaiu  and  the 

i-vat  ion  point  with  respeit  to  the  ^;h>ldow  bi'ondar  i es . 
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K 


ikp  jCos((u-  I t>j  I ) 


Jiji 


f E-wave'^ 

for  S ')  . 

. H-wave 
V .y 


(70) 


whoro  X is  defined  in  (12)  and  path  I'j  is  shown  in  Ei^-  12a.  In  order  to 
construct  (70),  \p  is  replaced  by  w in  (11)  for  reasons  which  will  become 
clear  later. 

Next,  we  consider  u|  as  the  incident  field  impinging  on  the  top 
half-plane  2.  The  total  field  diffracted  by  this  half-plane,  i.e.,  u^, 
may  be  split  into 


C t S 

U2  = I'l  + U2  , (71) 

R t 

where  u^  is  the  scattered  field.  Since  u^,  given  in  (69),  consists  of  two 

s s 

parts,  we  may  use  the  superposition  argument  and  designate  u^^  and  as 

is  s 

the  scattered  fields  due  to  u and  Uj,  respec  t iv'e  1 y . u^  can  then  be 
written  as 


a _ s s 

“2  "22  “21 


(72) 


t-  t s 

The  task  is  to  determine  u^,  which  now  consists  of  three  parts  as  u^,  u'.^2 

s t.  s 

and  n2j.  Fields  Uj  and  u‘22  can  be  constructed  from  expressions  given  in 
Section  2 as  they  are  the  total  and  scattered  fields,  respectively,  from 
the  lia 1 f-p lanes  1 and  2 illuminated  by  the  incident  plane  wave  (68). 

These  fields  take  uniform  forms  similar  to  (18)  and  their  asvmptotic 
expressions  can  readily  be  derived. 

In  determining  '‘21’  compare  (70)  with  (21)  and  employ  equations 


similar  to  (40)  and  (41),  arriving  at 


Figure  I2.i  imi  b.  (a)  Contour  F for  (70)  and  steepest  descent  path  SDl’  for  (77). 
(b)  Contour  F..,  for  (7  3). 
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whert>  'ti|  is  replaced  by  il  ■ 0 in  (70)  and  P2  ‘ind  are  defined  in  Fig.  12. 
Comparison  of  (7i)  with  (41)  reveals  that  x(f-.'u)  can  be  interpreted  as 
the  pattern  function  of  a fictitious  line  source  located  at  the  edge  of 
half-plane  1.  In  contr.ist  to  P('jj),  which  was  assumed  to  be  regular  in  (41), 
x(f2.,ai)  is  singular  at  w = *11^  and,  therefore,  special  attention  must  be 
devoted  in  the  asymptotic  evaluation  of  (71).  Attempts  were  made  to  first 
approximate  the  inner  ig-integral  in  terms  of  its  asymptotic  expansion  and 
then  to  evaluate  the  w- J nt egrat ion . Due  to  the  nonuniform  nature  of  tiiis 
procedure,  proper  results  could  not  be  obtained  for  all  the  cases.  Since 
the  inner  , -integral  can  be  performed  uniformly  in  terms  of  the  Fresnel 
integral,  we  rmty  now  express  (73)  as 


s 


T,  + tT.  + T.,  + tT, 
1 2 14 


(74) 


where 


2 


'3  ! 


I Xj(flj.ii))|  ikp,cos(w+4)  ) 

+ ^ V ikscos(.3-il)^  2 2 dw  , 

F]  |^x^(ni.u,)J 


(75) 


and  furthermore,  from  (13)  and  (20),  we  have 


fl  U) 

Xj(Oj.ui)  » +CSC 2 


to  + <)) 

= +/2k()2  sin  — 2 

r 


(76a) 


(7bb) 


i; I JcJ/  ■ 


It  is  not  iceti  th.it  on  p;ith  fj,  x^.  'S’  we  1 1 -boh.ived  whereas  x.  possesses  a 
pole  sinpiilaritv  at  u>  = :i..  Therefore,  asvmptotie  evaluation  of  T,  and  T, 
follows  the  s.ime  steps  used  to  evaluate  (42b).  In  the  following  sections, 
we  devote  our  .ittention  to  the  asymptotic  evaluation  of  Tj  and  T.,,  which, 
in  conclusion,  .allows  one  to  determine  u^  from  (71). 

4.2  Tot.al  Field  when  li  . 

_ _ 1 2 * 

For  this  case,  where  the  pole  ui  = -J.  and  the  saddle  point  w = fi  are 
separated,  we  may  di-form  the  p.ith  1'^  into  the  saddle  path  SDP,  shown  in 
Fig.  12a.  taking  into  account  the  residue  contribution  of  the  pole  lo  = fi . , 
and  arrive  at 


1 ° 4ti 


'.V.g  ,) 

+ Aj('  j>w)  e e F(+-,.)  dm 

SDP  r 


•‘i  ^ ^2''. 

+ F ►2k(;^  sin • 


ikscos(ft.-..)  ikp„cos(r’.+y.,) 
e ' e ^ . (77) 


Since  at  m = the  relation  f, . ^ 0 holds,  one  can  then  employ  the  results 

r 

of  Sec.  3.2.1  and  determine  the  asymptotic  value  of  (77).  The  same  pro- 
cedure c.an  also  be  used  to  evaluate  T2  and  T.^  which  completes  the  derivation 
of  defined  in  (74).  l.'sinp,  the  procedure  described  in  connection 

with  (71)  and  (72),  one  can  also  determine  the  asymptotic  values  of  u^ 
and  The  final  result  for  the  total  field  u^  is  finally  giver,  as 


U2  = x(ri^,fOR(ks)x(il,'f2)R(kP2^  ^ 

where  x and  g(’)  are  defined  in  (12)  and  (48),  respectively.  Equation  (78) 
agrees  completely  with  the  result  obtained  using  GTD  formulation; 
furthermore,  and  x(f2,<(<2)  "lay  be  identified  as  Keller's  diffraction 


coefficients  at  the  edges  of  half-planes  1 and  2,  respectively. 
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4.J  Total  Kiold  when  and  = SI  (I'iS’.*  H'b) 

This  is  an  interest  inr.  case,  hecanse  from  CTI)  consideration  the  obser- 
vation anple  coincides  with  tlie  shadow  boundary  of  the  diffracted  rays 
or  it>inat  inj;  from  halt-plane  2.  In  this  case,  as  in  the  previous  case, 
the  pole  u)  = S2j  and  t lie  saddle  point  oi  = SI  are  separated.  One  c.in  then 
replace  Tj  and  T2  from  (75)  with  their  counterparts  from  (77).  Since  at 
(0  = T the  relations  f . = 0 and  E,  ^0  hold,  evaluation  of  T,  and  T,  from 
(77)  and  (75),  respectively,  requires  special  care.  This  is  done  hv  using 
the  procedure  described  in  Sec.  3.2.2.  Determination  of  'l^  and  T^  from 

(77)  and  (75),  respectively,  may  he  completed  using  the  analysis  of 

t s 

Sec.  3.2.1.  Finally,  Uj  and  u ^ can  be  constructed  from  the  solution  of 
the  half-plane  illuminated  by  a plane  wave  as  discussed  in  Sec.  2. 
l.valuating  all  these  components  and  adding  them  up,  we  finally  arrive  at 
the  following  total  field 

t 1 1 ' 

u^  =4  x(S7.,i2)g(ks  + kp2)  + ! Tx(fipn)x^(Si,0)  - ~ x’(n.,S7) 

L.  ^ 

• g(ks)g(kp2)  + 0(k  ^ ^)  (79) 


where  = T x(f7,,oj) 

1 r)UJ  1 


a 


and  that  Comparing  the  term  containing 


x'  with  the  one  containing  P'  in  (59),  we  may  once  more  interpret  x ‘is  the 
pattern  function  of  a fictitious  line  source  located  on  the  edge  of  half- 
plane 1.  expression  (79)  can  also  be  obtained  using  UAT  or  MSD  formula- 
tions in  two  successive  steps  in  a mechanical  fashion  from  (66)  and  (65), 
respec  t i ve 1 y . 


4.4  Total  Field  when  fl^  = Q and  Q '<  (FlS-  10c) 


This  is  another  Important  case  for  which  the  edge  of  half-plane  1 
coincides  with  the  direction  of  the  incident  shadow  boundary  of  the  inci- 
dent field.  In  contrast  to  the  last  two  previous  cases,  the  pole  w = 


izs: 
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and  tlu'  saddle  point  u = SI  eoiiK  ide  tor  the  situation  under  investigation 
in  this  section.  Hence,  vhen  deforming  into  SO?  for  Tj  and  T2,  defined 
in  (75),  the  integration  passes  into  a pr i nc ipa 1 -va I ue  integral  along  SUP 
plus  half  the  residue  contribution,  namely. 


ikp  cos(  ,.>+i}'  ) 

T,  = ,"<-s<-o.s(-Sl)^.  2 2 ^ 

I /.  ,T  I ''  1 1 

SUP  , , r 


1 411 

2 


1 ..  /. 


+ Y 2kg.,  sin  - — I e 


n + I2)  ikp2fos(<c(-'!.2^ 


(80) 


Although  at  the  saddle  point  where  to  = SI  the  relation  f, . 0 holds,  the 

r 

procedure  devehiped  in  Sec.  3.2.2  can  not  be  used  for  (80)  as  x is 
singular,  in  contrast  to  P which  was  well-behaved  in  (56). 

The  integrand  of  (80)  is  expanded  in  terms  of  l.aurent  series  as 


ikp  cos(u.+«i^^)  _ :Li  _ 

+ X.(Si,u)  e “ F(+f.)  = — ~ A.  + B.  + 0(d)  - U) 

1 1 uj  - si  1 1 

r r r 


(81a) 


where 


.. 

ikp2  1 

A.  = 2e  e '^F(+C.) 
r r 

r_  2 

ikp  I ^^i  _ 

'^i  = 2e  A , 

L . 


u)=Si 

1 


u)=Si 


(81b) 


(81c) 


Substituting  (81)  into  (80)  and  using  saddle-point  integration,  we  arrive  at 

ikp.,cos(Si-(J)„) 


T = +B^g(ks)  - 2 X - N 

r r L ^ 

+ o(k“^/^) 


^2'  Iks 
e 


(82) 


where  takes  the  following  asymptotic  expansion 
r 


3 -3/'’ 

Bj  =•  +2g(kp2)  x^(“.'l'2)  + 

r ' r I u)«Sl 
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Having  ck-U'rminoil  Tj  and  T,,  oiu*  can  tlicn  evaliiaLc  T,^,  and 
ii'j'j  in  the  same  manner  described  earlier  to  finally  arrive  at 


and 


t 1 , , ,,  , iks 

u , = - x(  ,‘{’.,)g(kp^)  e + 


[TXj.cn, + 2x' (fil,^>2)  I g(ks)g(kp2) 

-3/ 


+ 0(k 


(K3) 


where  x'i  ,^2^  Yu\ 


Again,  the  appearance  of  x'  is  ‘in  important 


U)=fi 


observation.  It  must  be  realized  that  for  the  case  under  study  in  this 
section,  the  fictitious  line  source  interpretation  does  not  apply  anymore, 
because  otherwise  the  term  x'  would  not  have  appeared  in  the  solution.  The 
complexity  stems  from  the  fact  that  the  diffracted  field  from  half-plane  1 
in  its  shadow  boundary  direction  is  not  a local  plane  wave. 

4.5  Total  Field  when  11.  = Si  = iji^  (Fig.  1 Od ) 

This  is  a unitiue  situation  as  the  incident  direction,  the  observation 

direction,  and  the  line  connecting  the  edge  of  the  half-planes  are  all 

aligned  with  each  other.  Since  the  pole  to  = coincides  with  the  saddle 

point  0)  = Si,  one  can  transform  Tj^  and  T^  from  (75)  into  the  form  presented 

in  (80).  At  ui  = Si.  the  relation  ^ ^0  holds;  therefore,  evaluation  of  T„ 

1 r 2 

follows  the  same  procedure  used  in  Sec.  4.4  for  evaluation  of  T^  from  (80), 

i . e . , 


T2  = - j Xj.(«,fi)R(ltP2)  + 2x;(Si,Sl)g(kP2)g(ks)  + 0(k’^^^) 


(84) 


Similarly,  evaluation  of  T.^,  given  in  (75),  requires  formulation  of  the 
same  step.s  employed  in  Sec.  3.2.2,  i.e.. 


T„ 


- 1 ^^2 

-y  TXj.(Si,Sl)g[kC‘i  +P2)1  " TX^(0,Si)g(ks)g(kp2) 


(85) 
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T^,  dt'linoii  in  (75),  may  bt>  obtaineci  usin^;  the  construction  of  Sec.  3.2.1 
to  obtain 


= Tx^(f),Si)j-(ks)x^(i7,SOg(kp2)  . (86) 

Since  u|  and  defined  in  (71)  and  (72),  respectively,  are  the  solutions 

due  to  plane-wave  illumination  of  half-planes  1 and  2,  respectively,  one 
may  use  (18)  to  arrive  at 


i k ( s+p  ) 

+ 

ik (s+p^) 


TXj.  (n,fi)g  I k(s  + p ) 1 

1 1/  O 

+ Tx^dl.fl)  e ‘ g(kp2) 


(87) 

(88) 


To  complete  our  evaluation  of  u^,  the  only  step  remaining  is  to  deter- 
mine Tj  from  (80).  Substituting  <ij^  = S2  into  (80),  one  may  express  T.  as 


411 


n,  , n\  lk(s+p  )cos((ii-Q) 

— 2 — l|/2kP2  .sin  — — — e dm 


SOP 


^ ik(s-l-p2) 
4 " 


(89) 


Introducing  the  change  of  variable 


t = /2  e^"/^ 


sin 


(90) 


into  (89),  one  obtains 


j ik(s+p^)  . ik(s+P2) 

T,  = -r  e - T~  <■' 

1 4 2n 


(|(t)F|c.-'"'Vk72  t) 


(91a) 


where 


Q(t)  - t"'(l  + it^/2)  ' = t“'  - t + O(t^) 


(91h) 
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Substituting  (‘ilb)  into  ('Jla),  one  tiotes  tiiat  in  addition  to  the  integration 
of  type  (58a)  with  n = 1,  one  must  also  evaluate  the  following  integral, 
viz., 


C ,,  -iir/4/, — \ -k(s+p)t“  -1  , . 

-t  f e /ko2  t ; e t dt  = I (s) 


(92) 


Differentiation  of  l(s)  with  respect  to  s yields 


I'  (s)  = -k 


-k(s+P2)t 

Pt  t|  e t dt 


- (s  + p^) 


(93) 


for  which  the  result  given  in  (58a)  for  n = 1 is  incorporated.  Since 
I(s)  ••  0 as  s ► “>,  one  finally  obtains 


1 ( s ) = /p  2 


do 


s /o  (o  + P2) 


i tan 


■1  r2 


(94) 


T|  may  at  last  be  evaluated  to  yield 


11  -1  f^2\ 

^1  = 4 2T,  /TI 


R(ks)g(kp  ) . (95) 

2 (s  + P2 ) 2 


I 

I 


Now  we  have  enough  information  to  construct  U2  from  (71),  (72)  and 
(74),  thus  arriving  at 

, ik(s+p„)[ 


1I2  “ e 


2^'i  . 1 -1  ri 


iks 


^ + Ylj-  tan  y—  + 1 g(ks  + kp2>  + e ‘ g(kp2>  1 

i + g(ks)R(kp)l-  27T+ pT)' 

+ 0(k"’^^)  , 


(96) 


J 
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will*  rt' 


X (■  ,, 

r sin  Q 


) = 


cos  ii_ 

9 

2 sin-  U 


(9  7) 


This  is  an  important  result  and  can  not  be  obtained  via  simple  application 
of  rav  techniques. 

4.6  Compa r 1 son  wi  th  Other  Tec iin  iques 

In  this  section,  we  apply  the  formulation  of  uniform  techniques, 
discussctl  in  Sec.  i.  i.l,  to  the  geometry  of  staggered  parallel  plates. 

Results  obtained  by  a cascading  approach,  in  which  the  fields  scattered 
by  the  first  half-plane  are  directly  used  as  the  incident  fields  for  the 
second  half-plane,  are  compared  with  those  given  in  the  previous  sections 
using  spectral  analysis. 

For  the  situation  shown  in  Klg.  10a,  where  both  the  edge  of  half-plane 
2 and  the  observation  point  are  outside  the  transition  regions,  GTD,  UTD, 

MSD  and  L'AT  give  the  same  asymptotic  result,  and  are  in  complete  agreement 
with  (78).  This  is  the  simplest  situation  for  which  even  the  GTU  formulation 
provides  an  adequate  result. 

In  the  case  shown  in  Fig.  10b,  the  edge  of  half-plane  2 is  outside 
the  transition  region  and  the  observation  point  is  on  the  shadow  boundary 
of  half-plane  2.  The  GTD  formulation  can  be  used  to  determine  the  field 
diffracted  from  half-plane  1,  but  this  construction  fails  to  provide  the 
correct  result  for  the  field  diffracted  off  half-plane  2 along  the  shadow 
boundary  direction.  It  actually  gives  an  infinite  field  which  Is 
physically  unrealizable.  At  this  point  one  might  be  tempted  to  employ 
UTD  formulation,  as  in  (64),  since  it  would  predict  a bounded  (finite) 
result.  It  can  be  easily  verified  that  UTD  formulation  does  not  give  the 
correct  asymptotic  result  when  compared  with  the  exact  asymptotic  solution 
given  in  (97).  In  fact,  UTD  formulation  fails  to  predict  the  term  x'  Iti  (97). 


'4. 
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Two  success  i VO  app  1 i cat  ii>ns  of  MS!)  and  IIAT  formulations  do  indood  provide 
the  correct  result  .ind  demon;;trate  perfect  agreement  with  (97). 

For  tlie  situation  shown  in  Fij>.  lOc.  wliere  the  ed^i’  of  half-plane  2 
lies  at  the  shadow  boundary  of  half-platie  1 and  the  observation  point  is 
positioned  away  from  the  transitimi  region,  the  GTl)  formulation  obviously 
tails  to  determine  the  field  correctly.  it  is  then  logical  to  consider 
whetlier  or  not  L'TD,  MSI)  or  CAT  can  be  applied  twice  in  succession  in  a 
cascading  fashion  to  derive  the  final  result.  Since  the  incident  field 
on  the  half-plane  1 is  a plane  wave,  application  of  each  of  the  three 
formulations  gives  Sommerfeld's  solution  for  the  field  diffracted  bv 
half-plane  1.  We  can  readily  show  that  this  field  Lakes  the  following 
form  at  the  edge  of  half-plane  2 

Oj  = j + ry^(il,fl)g(ks)  + . (98) 

Field  (98)  can  now  be  viewed  .as  the  incident  field  on  half-plane  2. 
Hence,  straightforward  (mechanic.il)  application  of  UTD,  MSI)  and  UAT  would 
hopefully  provide  the  field  diffracted  by  half-plane  2 at  the  observation 
point.  After  some  s imii  1 i f icat i on , this  field  takes  the  following 
asymptotic  form 

t 1 i ks  - 1 / a 

u,.  = -r  y(U,ifi  ) c ■ g(kp  ) + t x (il,fl) g(ks ) x (SI,  !>., ) g (kp . ) + 0(k  “■) 

^ (99) 

Comparison  of  (99)  with  the  correct  asymptotic  expansion  given  in  (83) 
reveals  that  (99)  is  not  complete  and  misses  the  term  corresponding  to 
2x  ' (il,  02  ^ other  words,  straightforward  (mechanic.al  ) application  of 

UTD,  MSI),  and  UAT  does  not  provide  the  correct  asymptotic  result  when  the 
edge  of  half-plane  2 lies  at  the  shadow  boundary  (transition  region) 
of  half-plane  1.  This  important  observation  was  first  m.ide  by 
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Jones  1 2]  1 for  the  liAT  f ormii  1 .i  t i on . fiiTo  wo  h.ivi'  come  to  a similar 
important  observation  that  the  meih.inical  application  of  UTl)  and  MSD 
does  not  provide  the  corret  t result  either.  Fioersma  [ 22|  and  Lee  and 
Boersma  [ 19l  resolved  tin.’  difficulty  related  to  I'AT  by  usin;;  the  argument 
that  the  diffracted  field  (98)  from  half-plane  1,  which  varies  rapidly  in 
the  transition  region,  does  not  comply  with  a ray  field  behavior.  Hence, 
this  field  cannot  he  used  directly  (mechanically)  as  the  incident  field 
via  the  application  of  the  I'A'I'  formulation.  They  expanded  (98)  in 
terms  (if  an  infinite  summation  of  cylindrical  and  plane  waves  and  applied 
the  HAT  formulation  to  each  term  of  expansion  separately.  Their  final 
uniform  result  agrees  with  Jones'  uniform  result  and  with  that  given  in  (83). 

The  difficulties  with  the  straightforward  (mechanical)  application  of 
LTD,  MSD  and  HAT,  .as  expressed  in  the  preceding  paragraph,  also  exist  for 
the  case  shown  in  Kig.  lUd.  In  this  case,  both  the  edge  of  half-plane  1 
and  the  observation  point  lie  at  the  shadow  boundaries.  It  is  straight- 
forward to  show  that  the  "mechanical"  application  of  the  uniform  theories 
gives 


'*2  " } + 0(k  *)  (100) 

for  the  dominant  term  of  the  total  field.  However,  the  correct  dominant 
behavior  is  given  first  in  the  r.h.s.  of  (9h) , whicFi  is  obviously 
different  from  (100).  As  a m-atter  of  fact,  as  p (9h)  gives 

(1/2)  exp(iks  + ikp)  for  the  domin.ant  term,  i.e.,  50  percent  different 
from  the  "mechanical"  result,  (100),  of  the  uniform  theories.  Boersma  [22] 
and  Lee  and  Boersma  [19]  were  able  to  show  that,  bv  carefully  analyzing 
the  problem  in  the  m.anner  described  earlier,  they  could  obtain  the  complete 
result  agreeing  with  (96).  Howevt*r,  tlie  same  analysis  has  not  yet  been 
performed  for  the  MSD  formulation. 
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5.  Sl'KCTRjM.  DOMAIN  API'ROACH  TO  VERIKICATION  AND  RKF I NKMENT 
OF  ASYMPTOTIC  SODUTIONS 

).l  Int  roclm  L ion 

One  ol  till'  most  challenging  problems  in  the  solution  ol  high- 
frequency  scattering  analyses  is  the  establishment  of  the  acenrarv  of 
the  results  and  the  refinement  of  the  solution  when  the  need  for  its 
improvement  is  clearly  indicated.  The  difficulty  in  verifying  whether 
the  asymptotic  expression,  typically  derived  from  the  rav  ipproaih,  does 
indeed  solve  the  boundary  value  problem  under  ci>nsidi-rat  ion  stem?  ,‘ritr.irilv 
from  the  fact  that  there  is  no  obvious  way  to  "build  in'  i hi  hi  o-  r 
conditions  in  solution  procedures  based  on  ray  methods.  \n‘  t ^ ' ' 
is  that  the  high-frequency  solutions  are  often  constnu  ti  1 ; ’ • • < .. 

far  fields,  whereas  the  application  of  the  boundary  conditi  O'-  * r‘ 

requires  the  near-field  information.  In  contrast,  the  integral  enuat i >n 
formulation  for  the  scattering  problem  is  based  direct  Iv  on  the  .i|.  li.ation 
of  the  boundary  condition  and,  consequently,  the  boundary  condition  ihecl 
is  redundant  for  this  approach.  However,  the  conventional  moment  method 
solution  of  integral  equations  is  limited  strictly  to  the  low  frequency 
and  resonance  regions  as  the  matrix  size  becomes  unmanageably  large  beyond 
the  resonance  region. 

In  this  section  we  will  briefly  outline  a spectral  domain  method  for 
bridging  the  two  approaches,  viz.,  the  Integral  equation  and  asymptotic 
techniques.  The  hybrid  method  has  the  desirable  feature  that  it  not  only 
verifies  the  accuracy  of  the  ray  solutions  but  provides  a systematic 
means  for  improving  the  solution  for  a large  class  of  problems  of 
practical  interest.  This  fact  will  be  illustrated  via  two  typical 
examples  given  in  this  section.  Other  cases  have  also  been  treated  and 
may  be  found  in  \2'i\  and  \2U\. 
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6.2  Do  VO  1 opmont  of  Spi.H'tr,il  Donioin  rormii  ion  i)  f t lie  Integral  Kciuation 
and  Its  Appro.x  ima  to  So  hit  ion 

The  key  to  corahining  the  asymptotic  solution  with  the  integral 
equation  formulation  lies  in  recognizing  the  tact  that  the  Fourier  transform 
of  the  induced  current  on  a scatterer  is  directly  proportional  to  the 
scattered  far  fiehi  and  that  a good  approximation  to  this  scattered  field 
is  often  availahle  from  any  number  of  asymptotic  methods,  e.g.,  nri)  or 
the  spectral  approach  discussed  in  preceding  sections.  To  take  advantage 
of  these  facts  we  choose  to  work  with  the  " Four ier- t ransf ormed"  or 
"spectral  domain"  version  of  the  integral  equation  r.ather  than  with  the 
I'onvent  iona 1 spatial  domain  counterpart  of  the  same  equation.  We  begin, 
however,  with  the  conventional  e 1 ec t r i c- f i e 1 d integral  equation  (K-equation) 
for  a perfectly  conducting  .scatterer: 

(T;  * j)^  = , C = Creen's  Dyadic  (101) 

where  J(r')  is  the  unknown  induced  surface  current  density,  the  subscript 
t signifies  the  tangential  component  of  the  field  on  the  surface  .S  of  the 
scatteriT,  is  the  incident  electric  field  on  the  scatterer,  and 
symbolizes  the  convolution  operation. 

As  a preamble  to  Fourier  transforming  (101),  we  first  extend  it 
over  all  space.  To  this  end  we  define  a truncation  operator  0(A). 


0(A) 


f A^6(r 


r^)  dr 


r E S 


(102) 


where  6 is  the  Dirac  delta  function.  Let  0(A)  bt'  defined  as  the  complemen- 
tary operator 


0(A)  = A - 0(A)  . 


(103) 


6] 


Then  (101)  ean  ho  n-wriltt-n  as 

0 * j = 0(-k‘)  + (G  * J) 


(104) 


for  .ill  spaco.  I'siiij'  the  dofiniLion  of  Lho  Kourioi  transform  introduced  in 
Soc.  2,  now  oxLondi'd  to  tho  j’.oiural  three-dimensional  case,  we  can  write 
(104)  as 


''■1  = -Pi  ♦ K nos) 

whore  F = F[  (G  * ( J)|  and  Ej  is  the  tr.insform  of  tlie  lanKc’nti.il  component 
of  the  incident  field  trunc.ited  on  S,  with  all  L r.ins f ormed  quantities  being 
denoted  iii  this  Section  bv  (>n  top.  N'ote  th.it  the  convolution  oper.itor  in 
the  integral  equ.it  ion  is  t r.insfornii-d  into  an  .ilgebr.iic  product  upon  Fourier 
transformation.  Note  also  th.it  (10'))  h.is  two  unknowns,  viz.,  .1  .ind  F,  which 
must  be  solved  for  s i mu  1 t .ineous  1 y . 

At  this  point  one  can  construct  the  solution  ot  (105)  in  at  least 
two  w.'iy.s.  The  first  of  these  is  based  on  an  iter.itive  procedure  th.at 
begins  with  an  approximation  to  ,1,  the  trtinsform  of  the  induced  surface 
curr«.'nt  , di'rived  from  the  .ipplication  of  somi'  asymptotic  procedure  to  the 
scattering  problem.  The  following  equation  is  t lu-n  employed  to  generate 
the  next  order  solution  and  the  procedure  is  repeated  until  convergence  is 
.ichlevi'd: 

j(n+l)  ^ + f.  r'|GJ^"^|  - 0!f“' |G,/"^  I (106a) 

and  in  the  space  domain 


(106b) 
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It  may  bo  shown  that  tho  combinatiun  of  t Ito  second  and  third  terms  in 
(106a)  represents  an  approximation  to  T derived  from  tlie  n''*^  approximation 
ot  J,  i.o.,  Also,  it  should  be  realized  tlint  the  inverse  operator 

G * is  algebraic  since  G itself  has  the  same  nature.  The  check  for  the 
boundarv  ci'nditiun  may  be  readily  applied  in  this  procedure  since 
■1  *[Gd^"^l',  the  last  term  inside  the  square  brai^ket  in  (106a),  represents 
the  tangential  T.-field  produced  on  the  surface  of  the  object  by  the 
indui-ed  current.  Ideal Iv,  this  field  should  equal  the  negative  of  the 
tangential  component  of  the  incident  E-field  E* . A comparison  between 
these  two  fields  immediately  reveals  the  extent  to  which  the  boundary 
condition  is  satisfied  on  the  surface  of  the  scatteri'r. 

A Second  approach  to  handling  (106)  would  u-  ti>  employ  the  Galerkin 
procedure  in  the  transform  domain.  One  may  write 

J + )■  C .1  (107) 

where  is  the  approximate  solution  derived  from  a suitable  asymptotic 

formula  for  the  scattered  field  and  .1^  represents  a set  of  basis  functions 

in  the  transform  domain.  Typically,  there  are  certain  angular  regions  in 

the  far  field  where  the  asymptotic  solutions  require  refinement.  One  may 

choose  to  concentrate  the  basis  functions  in  these  regions  in  the  transform 

domain.  Alternatively,  the  .1  's  could  be  chosen  as  the  transforms  of  a 

P 

suitable  set  of  basis  functions  in  the  space  domain,  and  the  location 
(support)  of  these  subdomain  basis  functioiis  m.iy  be  selected  to  coincide 
with  transition  regions  or  corners,  etc.,  where  the  canonical  solution  of 
the  asymptotic  solution  may  require  refinement. 

In  either  case,  the  problem  ot  determining  .1  may  be  reduced  to  that 
of  finding  the  unknown  coefficients  such  th.at  (107)  satisfies  (105). 
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The  C.'.ilerkin  proredure  provides  ;i  wav  for  acoomp  I i sii  i np  this,  as  wo  will 
soiMi  see.  This  toi'hnitpK'  also  has  the  adv.inta^o  that  tlio  otlior  unknown 
in  (103),  vi/.,  K,  is  ronvcn  i out  1 y oliminatod  from  this  fciu.ition  upon 
appl  icat ion  ot  Galorkin's  method.  We  demonstrate  this  fact  in  tlie  mani- 
pulat  ions  presented  beli'w. 

Suiist  i t ut  i ng  (107)  in  (105)  and  takinp,  a scalar  product  of  the 
resulting  equation  with  a set  of  suitable  series  of  testing  functions  W , 

q 

we  arrive  at 


y C "W  ,0,T  > = -<W  ,£,  ■ + -W  ,F- 
^ p q p q I q 


(108) 


where  "■ , ''  is  the  scalar  inner  product.  If  we  now  choose  W to  be  trans- 

q 

forms  of  functions  which  are  nonzero  only  on  the  surface  of  the  scatterer, 
then  the  scalar  product  <W  , F>  can  be  shown  to  vanish.  To  show  this,  one 

q 

uses  I’arseval's  theorem  and  transforms  the  scalar  product  of  W and  F in 

q 

terms  of  a similar  product  of  their  counterparts  in  the  space  domain. 

Since  the  inverse  transforms  of  and  F exist  in  complementary  regions, 
viz.,  on  the  surface  of  the  scatterer  and  in  tlie  region  complementary  to 
this  surface,  respectively,  one  finds  that  their  scalar  product  is 
identically  zero.  One  can  now  proceed  in  the  usual  manner  to  solve  for 
the  coefficients  by  solving  the  matrix  equation  represented  by  (108) 

with  the  second  term  in  the  r.h.s.  deleted.  It  is  evident  that  the  use 
of  this  method  would  he  practical  only  when  relatively  few  terms  are 
needed  in  (107)  to  modify  the  available  asvmptotic  solution;  however, 
this  is  typically  the  situation  for  many  problems.  It  should  also  he 
noted  that  (108)  represents  a dirc'ct  I'heck  on  the  satisfaction  of  the 
boundary  condition,  in  the  sense  of  moments.  The  choice  of  is 

governed  hy  the  locations  on  the  surfaci's  of  the  scatterer  where  these 


- — »,;•>  ••  .. 
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bouiulary  conci  i t i mis  aro  applied.  Tvpirally  those  will  be  the  zones  where 
the  asvmptot  io  solution  raif’ht  be  inacourate,  e.g.,  the  transition  rep, ion 
between  the  lit  and  shadow  regions. 


5.  ) 1 1 1 list  rat  i ve  bxamii  les 

We  will  now  illustrate  both  the  iterative  and  Calerkin  procedures  by 
nu'.ins  ot  two  examples  of  high-f reqiiencv  scattering.  For  the  first  example 
we  choose  the  problem  of  diffraction  by  a rectangular  cylinder,  a structure 
with  multiple  wedges.  The  geometry  of  the  problem  is  shown  in  Fig.  13. 

We  note  th.it  the  wi-dges  1,  2,  and  3 and  faces  A and  B are  illuminated  by 
the  incident  wave  while  the  faces  C and  D and  the  wedge  4 are  in  the  shadow. 
Let  us  first  consider  the  heimsphere  above  the  face  A.  The  contribution 
to  the  scattered  field  in  this  region  comes  from  the  two  lit  surfaces  and 
rays  diffracted  by  the  wedges  1,  2,  and  3.  To  derive  a first  approximation 
to  the  asymptotic  expression  for  the  far  field  scattered  by  the  cylinder, 
we  may  use  Keller's  coefficients  for  the  diffraction  by  the  three  wedges. 
From  wedge  1 we  have,  for  instance. 


- i k (-acosif^+bs  in(f^)  ikp 


0' 


, I 2 . 2tt  ikacosif  -ikbsinf) 

4 |—  sin  j-  e e 


zl 


2ti 


</p  /2-nk  I cos  -j  - cos  ^ (<p  - 'p^) 

2 . 2n  ikacosif  -ikbsinf 

— sin  - e e 

cos  ^ - cos  ^ (if  + ly) 


(109) 


where  0 •'  4’q  0 < if  < 3n/2.  Similar  expressions  can  be  readily 

obtained  for  the  wedges  2 and  3.  Note  that  the  scattered  field  from  wedge  3 
is  restricted  in  the  angular  region  0 < if  •'  ir/2  and  n < if  < 2n  for 
0 ^ n/2.  Hence,  in  the  upper  hemisphere  only  one-half  of  the  angular 

range  Is  illuminated  by  the  scattered  field  from  wedge  3,  if  Keller's 


FiKure  13.  I’of  1-l‘f  .inglf  ol  i n<- i dtiu-c  ; .i-:  shown,  wnd^es  1,  2,  and  3 

are  illuminated  while  wedKi-  5 is  in  the  d.irk. 
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lurmiila  for  ;in  inliniti-  wodj>,o  is  used  to  dorivo  an  a|)[>  rox  i nia  t i on  to  tfiis 
soattorod  fii’ld.  Another  ratlier  important  and  common  feature  of  the 
Keller  expressions  for  the  weilt’e-d  i f f rac  t ed  fields  as  ttiven,  for  instance, 
by  (109),  is  that  tltese  formulas  predict  fictitious  infinities  in  the 
sc.ittered  fit‘lds  at  thi'  shadow  anil  reflection  boundaries.  One  could 
completely  eliminate  these  intinities  by  employiiiH  the  spectral  concepts 
and  deriving  t lu‘  scatti-red  lar  field  from  the  transform  of  an  approximation 
t the  indiiceil  surface  current,  comprising  the  physical  optics  and  fringe 
currents.  This  transform,  being  associated  with  a function  with  finite 
support,  is  .ilwavs  bounded  and  consequently  free  of  the  singular  behavior 
prisent  in  (109). 

As  .in  alternative  we  may  also  employ  one  of  the  available  uniform 
theories  (17|,  [14]  that  jirovide  smooth  and  bounded  transition  through  the 

ri'f  lection  .and  shadow  boundaries.  It  is  fortuitous,  however,  that  the 
.aggregate  contribution  of  the  infinities  from  the  individual  wedges  cancels 
out  exactly  when  their  contributions  are  superimposed.  This  occurs 
because  of  the  unique  symmetries  of  the  geometry  of  the  rectangular  cylinder 
under  consideration.  Hence,  no  special  care  is  required  in  this  example 
.at  the  transition  regions.  The  diffracted  far  fields  computed  by  using 
the  Keller  formulas  for  wedge  diffraction  are  shown  in  Fig.  14.  It  is 
evident  that  the  pattern  is  discontinuous  at  0°,  90°,  180°,  and  270°. 

This  behavior  is  attributable  to  the  use  of  infinite  wedge  diffraction 
coefficients  that  produce  a nonzero  scattered  field  only  in  the  region 
external  to  the  wedge  and,  hence,  produce  discontinuous  fields  supported 
by  induced  surface  currents  that  extiuid  to  infinity  along  the  wedge  surfaces. 

To  refine  this  far-field  asymptotic  si'lution,  we  m.ay  now  proceed  to 
.apply  the  iterative  procedure  developed  in  Sec.  5.2.  To  this  end,  we 
first  introduce  in  Fourier  transform  the  variable  it  = k.  cos  i|i  and  express 
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till'  t.ir  Si  .It  Lffi'il  lifid  in  t hr  horn  i splie  r i c' a 1 rogion  dofinod  by  0 ' n 

in  triTis  of  1.  Wo  omp  I I'v  .an.i  1 y t ira  1 oont  i nna  t i on  of  tho  oxpressiuns  for 
till'  wodgi'-d  i t I rar  t rd  lii'lds  to  dotormino  tho  Konrior  transform  in  liie 
ranpo  i|  • k bv  snbst  i t nl  i np  appropriate  complex  values  for  i- . One  of 
till'  I'liief  advatUajtos  of  using  the  available  expressions  for  the  asvmptotic 
solution  as  a starting  point  for  the  iterative  procedure  is  that  the 
approximate  analytical  I'xpression  is  convenient  for  estimating  the  scattered 
f.ir  field,  both  in  the  visible  and  invisible  ranges.  By  Fourier  inversion 
of  the  scattered  E-field  at  infinity,  we  can  derive  the  tangential  E-field 
on  a planar  surface  tangential  to  the  face  A.  If  this  computed  near  field 
were  to  satisfy  the  boundary  condition  that  the  tangential  E-field  on 
surface  A equal  the  negative  of  E^,  and  the  similar  situation  is  repeated 
for  other  faces,  we  would  conclude  that  the  solution  so  derived  is  an 
accurate  one.  However,  we  would  not  expect  that  to  be  true  for  the 
solution  represented  by  (109)  since,  as  pointed  out  earlier,  this 
expression  produces  glitches  at  some  angles  of  observation.  Nevertheless, 
this  solution  provides  a very  good  zero  order  estimate  for  which  is 

readily  derived  by  taking  the  transform  of  the  near  field  derived  in  the 
plane  of  the  surface  A,  after  deleting  the  portion  of  the  field  that 
corresponds  to  the  surface  A of  the  scatterer,  and  rejjeating  the  same 
procedure  for  the  other  faces  as  well.  In  following  this  procedure,  we 
effectively  compute 


is  Fourier  inverted  four  times  to  calcnl.ite  the  tangential  H-field  on  the 
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tour  f.u'os  of  the  cylinder  and  the  surface  current  on  these  faces,  obtained 
from  (106b),  is  Fouric'r  transformed  .igain  to  derive  the  far-field  pattern. 
The  iterated  far-fii'ld  pattern  .1^'^  is  shown  in  Kig.  13;  the  disappearance 
of  the  d isi-ont  inui  t ies  at  0,  ti/2,  '3n/2,  etc.,  is  immediately  evident  from 
this  plot.  This  result  has  also  been  verified  by  a few  other  workers  who 
have  followed  different  procedures  than  those  outlined  liere  [251.  Recall, 
howi'ver,  that  the  method  outlined  here  provides  a convenient  "built-in" 
check  for  the  satisfaction  of  the  boundary  condition  and  an  independent 
check  is  not  altogether  necessary  to  establish  the  accuracy  of  the  solution. 

The  next  example  we  discuss  is  that  of  a smooth  convex  surface  with 
no  wedges  a circular  cylinder.  One  of  the  important  attributes  of  this 
canonical  geometry  is  that  it  permits  convenient  comparison  with  the  exact 
series  solution  available  for  the  representation  of  scattered  fields  from 
this  structure.  The  geometry  of  the  problem  is  shown  in  Kig.  16.  We 
consider  the  case  of  an  E-polarized  wave  incident  from  = 180°. 

The  first  step  in  attacking  the  problem  is  to  use  a geometric.il 
optical  (60)  approach)  |2A)  to  derive  the  far-scattered  field.  When  this 
is  done,  one  obtains  the  dotted  curve  in  Kig.  17  which  also  exhibits  the 
exact  series  solution  for  the  scattered  field  as  a solid  curve.  It  is 
evident  that  in  the  range  -60°  < <\i  < 60°  the  CO  solution  is  not  adequate. 
This  is  not  totally  surprising  since  it  is  well-known  that  creeping-wave 
contributions  need  to  be  included  in  the  scattered  field  expression  in  the 
shadow  and  transition  regions.  Rather  than  following  this  procedure,  we 
will  now  show  how  the  Galerkin  method  can  be  readily  and  conveniently 
applied  to  this  problem  to  derive  an  accurate  solution. 


* Numerical  experiments  have  shown  that  very  good  f.ir-field  results  are 

obtained  even  when  one  skips  the  first  step  .and  uses  d^*'^  itself  to  derive 
the  surface  current  using  (106b)  and  d^^^  by  the  Fourier  transformation  of 
this  current. 
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Figure  16.  Diffraction  by  a circular  cylinder  illuminated  by 
E-wave  incident  along  the  x-axis. 


an 


7 3 

lo  tills  LMid  wf  consider,  as  a first  step,  the  beliavior  of  the 
■scattered  fitld  on  a surlaee  erected  in  juxtaposition  to  the  cylinder  at 
tin.'  iH'int  X = a,  the  farthest  point  awav  from  the  incident  field.  In  tfie 
deep  shadow  region,  say  jyl  < 2,  we  exjiect  the  sc.ittered  lield 
to  b(  a verv  gooil  approximation.  On  the  other  hand,  when  we  go  far  onto 
t lu'  lit  rcj’,  ii'ii  on  this  surface,  say  for  |y  h , we  expect  the  to  lie 
descrihed  adecpiately  hy  the  flO  formal. is.  It  we  had  a good  estimate  of 
the  scatteri'd  lield  behavior  in  the  transition  region  2 < |yl  •"  6,  we 
wi'uld  lu-  able  to  get  a good  representation  of  the  excess  scattered  field 
(over  and  abovi-  the  (It)  field)  on  the  entire  surface  at  x = a.  We  should 
then  be  ible  to  compute  the  field  radiated  in  the  r.h.s.  of  the  cylinder 
hv  this  i-xcess  field  using  the  concept  of  Huyghen's  source  ,ind  use  this 
radiated  lield  to  fill  in  the  gap  between  the  (10  pattern  and  the  true 
|iattern  shown  in  Fig.  18. 

To  derive  the  E'*’  field  in  the  transition  region,  we  first  interpolate 

7. 

till-  m.ign  i t ude  of  this  field  from  !'.  * ,at  |y|  = 2,  to  0 at  y = (i  and  t he 

[ihase  from  :r  at  |y|  = 2,  to  the  CO  phase  at  v[  = C . Next  we  introduce  a 

sit  of  basis  functions,  with  undetirmined  coef f 1 c ient s , to  describe  the 

correi  t ion  to  the  interpolated  1.'^  I ield  at  the  plane  x = a,  2 < |y|  < b. 

To  determine  th»-se  coe  I f i c i en  t s , we  apply  tlu-  concept  of  Calerkin's 

method  in  the  snectral  doin.iin  as  briefly  outlined  in  Sec.  5.2.  in  the 

i-xami>le  being  considered  here,  the  tr.insforms  of  the  basis  functions  play 

the  role  of  .1  in  (107),  and  tlu-  zero-order,  sc.attered  far  field  is 

P 

obtaiiu-il  by  adding  tlu-  contributions  of  CO  and  the  approximate  t-xcess 
E]^  fi(-ld  d(-rived  from  tlu-  interpolation  procedure  just  descrihed. 


The  choice  of  Llie  Lest  inf,  functions,  in  (108),  is  suSL^t^sted  bv  the 
fact  tli.it  tile  error  in  the  h igh- f requi'ncy  asymptotic  solution  is  mostly 
c'oncent  rateii  around  the  transition  region  on  tile  surface  of  the  cvlinder, 
i.e.,  in  the  neighborhood  of  the  junction  between  the  lit  and  shadow 
regions.  Thus,  a suitable  choice  for  the  testing  fumtions  would  he  to 
locati'  them  at  the  transition  region  as  shown  in  Fig.  18  where  the  location 
of  the  basis  functions  is  also  shown.  Note  that  we  need  not  be  restricted 
in  our  choice  for  the  location  of  these  functions  by  demanding  that  Lhc“v 
have  a common  support,  although  tliis  is  almost  always  the  case  in  the 
lonvent iona 1 moment  or  Galerkin  methods.  We  may  also  note  from  Fig.  18 
that  the  shape  of  the  basis  and  testing  functions  are  both  C.aussian. 

Since  we  are  dealing  with  transforms,  this  choice  is  not  only  convenic-nt 
for  deriving  the  Fourier  transformations  J and  W , but  is  also  desirable 

p Cl 

from  a numerical  point  of  view  because  the  transforms  are  not  oscillatory 
as  they  would  be  for  a pulse  or  triangular  basis.  This  feature  is  Important 
when  numi'rically  computing  the  scalar  products  , needed  for  the 

determination  of  the  unknown  coefficients  C . 

P 

Only  a few  unknowns  C ('3  to  7)  are  needed  to  derive  an  accurate 

P 

solutitjn  for  both  the  radiated  far  field  and  the  surface  current  on  the 
cylinder.  The  accuracy  itself  can  be  verified  by  computing  the  tangential 
E-field  on  the  surface  via  Fovirier  inversion  of  the  hemispherical  fir- 
field  pattern  centered  around  the  point  to  be  tested.  This  procedure  is 
also  usiui  to  compute  the  surface  current  distribution  from  the  knowledge 
ot  the  scattered  11-field  at  large  distances.  Of  course  an  independent 
check  Is  avaii.'bli-  for  this  problem  via  the  exact  series  solution.  A 
comparison  <>l  the  O.ilerkln  solution  and  the  exact  series  solution  is 
shown  in  Figs.  19  and  20  to  Illustrate  the  highly  accurate  nature  of  the 
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C.'ilorkin  solution.  in  the  solution  is  su  ri>r  i.s  inf>  I y accurate,  except 

for  a sliglit  erri>r  in  the  transition  region,  even  without  the  fialerkin 
refinement,  as  evidenced  by  the  dotted  curve  in  Fig.  20  which  exhibits  this 
case . 

The  met  liod  fust  described  couhi  easily  be  employed  for  more  complex, 
nonseparab 1 e shapes  as  long  as  the  smooth  convex  nature  of  the  surface  is 
preserved.  A gene r;i  1 i za t i on  to  even  more  complex  geometries,  e.g., 


combination  of  curved  surfaces  and  wedges, 


also  appears  possible. 
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